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at  4K  we  have  deduced  information  about  the  alloy  potential  fluctuations  in 
MBE  fabricated  grown  at  different  substrate  temperatures. 
A  comparative  study  of  RS,  transmission  electron  microscopy  (TEM)  and 
conductivity  has  been  performed  on  molecular  beam  deposited  microcrystallir.e 
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about  the  surface  e'lctric  fields  in  heavily  doped  <100>  n-GaAs  has  been 
obtained  by  stifding  the  symmetry  forbidden  polarization  configurations  in 
RS.  We  have  performed  a  detailed,  comprehensive  investigation  of 
polish-induced  surface  strain  in  <100>  and  <111>  GaAs  and  InP.  Considerable 
information,  has  been  obtained  about  surface  strain,  strain  skin-depth  and 
diffusion  of  damage  into  the  sample.  • 

A  photoreflectance  (PR)  study  of  Hg^_^Cd^Te  fabricated  by  different  • 
growth  conditions  has  demonstrated  the  utility  of  this  method  in  evaluating 
crystal  quality.  In  addition,  by  investigating  the  pump  chopping  frequency 
dependence  of  the  PR,  information  was  obtained  about  surface  (interface) 
traps.  The  interface  properties  of  Hg.  ,Cd-  ,Te/SiO  were  studied  by 
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photo-and  bias-induced  charging.  It  was  possible  to  detect  very  slow  ' 

interface  states  and  to  gain  information  about  their  energy  distribution.  ■, 
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I .  INTRODUCTIOIL 


Raman  scattering  (RS)  is  an  extremely  usetul  tool  for  studying  semicon¬ 
ductors  either" in  bulk  or  thin  film  form.  RS  depends  on  the  electron-phonon 
interaction  and  therefore  is  an  important  complement  to  other  spectroscopic 
techniques  (ellipsometry,  luminescence,  modulated  spectroscopy,  various 
photo-electron  spectroscopies,  etc.)  which  are  functions  of  the  electron 
states  of  the  solid.  RS  is  one  of  the  most  important  methods  for  studying 
the  lattice  vibratiion  of  a  semiconductor  and  interactions  with  other 
excitations  such  as  plasmons,  magnons  etc.  Since  lattice  vibrations  are 
very  sensitive  to  local  environments,  RS  can  yield  information  about  the 
semiconductor  on  the  scale  of  a  few  lattice  constants.  Thus  it  is  ideal  for 
investigating  the  nature  of  microstructural  geometries,  i.e.  structures 
which  are  of  order  10  -  1000  8  in  dimensions.  Also  since  RS  is  a 
second-order  process  it  contains  important  symmetry  information  not  avail¬ 
able  from  first-order  optical  interactions.  In  RS,  both  polarization 
selection  rules  as  well  as  peak  positions  are  sensitive  to  perturbations, 
both  internal  and  external,  such  as  strain,  electric  fields,  temperature, 
etc.  Furthermore  by  using  various  excitation  lines  with  different 
wavelenghts  it  is  possible  to  (a)  perform  depth  profiling  measurements  or 
(b)  excite  resonance  interactions. 

One  of  the  most  significant  aspects  of  some  of  our  initial  results  has 
been  the  demonstration  the  a  "spatial  correlation"  model  (SCM)  with  a 
Gaussian  correlation  function  is  a  very  general  approach  to  q-vector  relaxa¬ 
tion  in  semiconductors  due  to  a  number  of  microstructural  effects.  This 
model  has  already  been  used  to  account  for  finite  size  effects  in 
microcrystalline  geometries  but  the  most  recent  experiments  show  that  it  is 
also  applicable  to  ion-dairiage  material  and  alloy  semiconductors.  The 
detailed  emalysis  of  the  Raman  lineshape  yields  a  physically  significant 
parameter,  i.e.  the  "spatial  correlation  length  which  can  be  related  to  the 
size  of  the  microcrystalline  region,  the  size  of  the  ion-damage  region  in 
ion-damaged  (and  possible  annealed)  material  or  the  size  of  the 
compoistionally-ordered  section  of  alloy  semiconductors.  The  ability  to 
extract  this  microscopic  quantity  from  macroscopic  measurements  opens  up  new 
areas  of  investigation  that  are  of  significance  from  both  a  fundamental  and 
technological  point  of  view.  Using  the  SCM  we  have  analyzed  the  RS  of  the 
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alloy  AIq  grown  under  different  conditions  including  sub¬ 

strate  temperature.  Evidence  for  alloy  ordering  is  found  at  a  particular 
substrate  temperature.  Molecular  beam  deposited  microcrystalline  Si  and  Ge 
has  been  studied  by  RS,  Transmission  Electron  Microscopy  (TEM)  and  conduc¬ 
tivity.  It  is  found  that  RS  can  detect  point  defects  not  observed  by  TEM. 
Information  about  the  electric  fields  in  the  space-charge  region  of  heavily 
doped  <100>  n-GaAs  has  been  obtained  from  symmetry  forbidden  polarization 
selection  rules.  We  have  performed  the  firs  comprehensive  study  of 
polished- induced  strain  in  <100>  and  <111>  GaAs  and  InP.  A  detailed 
lineshape  analysis  of  the  RS  has  yielded  important  information  about  surface 
strain,  strain  skin-depth  and  the  diffusion  of  damage  into  the  material. 

In  the  second  phase  of  our  program  we  investigated  (a)  the  excitation 
wavelength  and  pump  chopping  frequency  dependence  of  photoreflectance  (PR) 
in  Hg^  ^Cd^Te  and  (b)  the  properties  of  the  Hg^  ^Cd^  ^Te/SiO^  -Photox 
interface  using  photo-and  bias-induced  charging.  The  former  experiment 
yielded  information  about  the  quality  of  Hg  _  Cd  Te  produced  by  various 
growth  conditions  as  well  as  surface  states.  The  latter  work  gave  signifi¬ 
cant  insights  into  the  natures  of  traps  at  the  Hg^_j^Cd^Te/Si02  -Photox 
interface. 
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II,  EXPERIMENTAL  RESULTS 

A.  Raman  Spectroscopy  Study  of  the  Effect  of  Growth  Conditions  on  Alloy 

Potential  Fluctuations  in  Al.  __As/InP 

0.48  0.52 

We  have  measured  300  K  Raman  scattering  (RS)  from  several  samples  of 

<100>  AIq  431‘^Q  grown  by  molecular  beam  epitaxy  with  substrate 

temperatures,  T  ,  in  the  range  495°C  <T  <550°C^.  The  4°K  photoluminescence 
s  s 

(PL)  from  these  samples  has  also  been  studied  .  The  details  of  the  longitu¬ 
dinal  optic  phonon  RS  line-shape  are  interpreted  in  terms  of  a  "Spatial 
Correlation  Model".  At  =  515°C  there  is  a  significant  decrease  in  both 
RS  and  PL  linewidths  indicating  that  under  this  growth  condition  there  is  a 
decrease  in  the  alloy  potential  fluctuations,  i.e.,  a  tendency  towards 
ordering.  The  "Spatial  Correlation  Model"  line-shape  analysis  is  compared 
with  a  recently  proposed  "Discrete-Continuum  Interaction  Model",  i.e., 
Fano-type  interaction. 

The  ternary  alloy  Al^In^_^As  is  important  for  both  fundamental  and 

technological  reasons.  With  regard  to  the  former  considerations  this 

material  should  exhibit  bond- length  (and  bond-angle)  distortions  similar  to 

those  found  recently  in  Ga^In^_^As^.  The  resultant  internal  strain  is  found 

to  be  larger  for  the  disordered  phase  in  relation  to  the  ordered  phase  and 

thus  may  serve  to  stabilize  the  ordered  phase^  Orderings  in  the  systems 
6  7  8  9 

GaAlAs_  ,  InGaAs  ,  Ga  AsSb  and  SiGe  have  recently  been  reported.  The 

^  ^  ^  IQ 

role  of  the  growth  temperature  has  been  discussed  by  Srivastava  et.  al  and 

Mbaye^^  et.  al.  On  the  basis  of  their  first-principles  total-energy  and 

phase-diagram  calculations  these  authors  have  proposed  that  whereas  at  high 

temperatures  {T>T^)  such  systems  are  stable  in  the  disordered  phase  at  lower 

temperatures  (T<T^)  an  ordered  phase  or  clustered  (pseudo-ordered)  phase  can 

become  the  ground  state.  Ichimura  and  Sasaki's  calculation  shows  that  there 

12 

is  a  preference  for  ordering  in  the  III- IV  ternary  alloys.  In  the  area  of 

technology,  devices  using  the  Al^  43^'^Q  53^^  heterostucture 

include  high-mobility  modulation-doped  structures  ,  double-heterostructure 

lasers  which  emit  at  1.55  um^^,  quantum  will  structures  which  show  promise 

13 

for  lasers  emitting  from  less  than  1.3-1.65  Min  »  high-speed  long-wavelength 

14 

quantum-well  avalanche  photodiodes  ,  etc. 
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We  have  studied  300K  Raman  scattering  (RS)  measurements  of  a  number  of 
samples  of  <100>  Al^  43!"^  grown  by  MBE  with  substrate  temperatures 

in  the  range  495°C  <  T  <  550°C.  It  should  be  noted  that  Al„  ,„In„  ^_As  is 
approximately  AlInAs^.  For  comparison  purposes  we  have  also  studied  the  4K 
photoluminescence  (PL)  of  these  same  samples.  It  has  recently  been  suggest¬ 
ed  that  and  RS^^  cam  be  used  to  study  the  microscopic  alloy 

potential  fluctuations  in  alloy  semiconductors.  There  is  a  strong  correla¬ 
tion  between  the  details  of  the  RS  line-shape  and  the  PL  linewidth.  At  T  = 
Q  s 

515  C  there  is  a  significant  decrease  in  both  the  RS  and  PL  linewidths, 
indicating  that  under  this  growth  condition  there  is  a  decrease  in  the  alloy 
potential  fluctuations,  i.e.,  a  tendency  towards  ordering. 

The  RS  of  this  alloy  system  exhibits  a  two-mode  behavior  corresponding 
to  "AlAs-like"  and  "InAs-like"  lattice  vibrations.  The  first  order  longitu¬ 
dinal  optic  (LO)  phonon  Raman  lines  exhibit  a  low-frequency  asymmetrically 
broadened  line-shape.  We  have  analyzed  this  line-shape  in  terms  of  the 
"Spatial  Correlation  Model"  (SCM)^^.  The  line-shapes  of  the  "AlAs-like" 
mode  for  all  samples  and  the  "InAs-like"  mode  for  T^  =  515°C  can  be  fit  by 
the  SCM  using  only  one  adjustable  parameter,  i.e.,  the  "Spatial  Correlation 

Length",  However,  for  the  "InAs-like  mode  for  T  f  515°C  both  I  and 

s 

the  intrinsic  linewidth  (r^)  must  be  varied  to  achieve  a  fit  (two  parameter 

fit).  This  indicates  that  for  these  growth  conditions  a  new  scattering 

mechanism  has  become  available  for  this  mode.  The  recent  suggestion  that  a 
21-22 

Fano  model  (two-parameter  fit),  rather  than  the  SCM,  can  be  used  to 

account  for  the  LO  phonon  line-shape  of  alloy  semiconductors  will  be  dis¬ 
cussed. 

The  details  of  the  sample  growth  procedures  are  given  in  Ref.  13.  The 
AlInAs  epilayers  were  typically  2500A.  The  RS  measurements  were  performed 
at  300K  using  the  5145A  line  of  an  Ar-ion  laser.  The  scattering  geometry 
was  x(y,z)x  where  x  =  <100>,  y  =  <010>  and  z  =  <001>.  For  this  configura¬ 
tion  only  LO  phonon  scattering  is  allowed  in  the  first-order  RS.  To  mini¬ 
mize  heating  effects  a  line  focus  ( lOOum  x  3  mm)  was  used  for  the  laser  beam 
(typical  power  of  1  watt)  and  the  sample  was  thermally  bonded  to  a  heat 
sink. 

Shown  in  Fig.  1  is  the  300°K  RS  of  Al^  52^^^^^^  sample  C  (T^  = 

515°C).  The  spectrum  is  typical  of  all  the  samples  used  in  this  study.  The 
large  feature  at  about  240  cm  ^  is  identified  as  the  "InAs-like"  LO  phonon 
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is  the 


mode  and  is  designated  LO^  while  the  peak  at  about  365  cm 

"AlAs-like"  LO  mode,  labelled  LO  .  The  frequencies  of  both  these  modes 

14 

agree  well  with  other  results  .  The  weak  broad  structure  around  180  cm 
is  probably  related  to  disorder-activated  zone-edge  modes . 

19 

In  the  SCM  the  Raman  intensity  1(a))  at  frequency  o)  is  given  by 

St 


If'-') 


I 


(1) 
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where  BZ  denotes  the  Brillouin  zone,  «,is  the  "Spatial  Correlation  Length" 
and  Fq  is  the  intrinsic  linewidth  of  the  relevant  end  point  material.  For 
the  dispersion  aj(q)  of  the  LO  phonon  we  take  the  analytical  relationship 
based  on  the  linear  chain  model: 


aj^(q)  =  A  +  {A^  -  B(  l-cos(nq)  ]  (2) 

4-2  9-4 

For  the  "AlAs-like"  modes  we  take  A  =  8.2x10  cm  and  B  =  2.23x  10  cm 

4  -2 

while  for  the  "InAs-like"  modes  the  parameters  are  A  =  2.83x10  cm  and 
8  -4 

B  =  3.4x10  cm 

Shown  by  the  dotted  lines  in  Fig.  2  are  the  experimental  data  for  the 

LO,  and  LO.  Ramein  bands  of  Al„  -oln,,  ^_As/InP  samples  A,  C  and  F.  Represen- 
i  ^  U . 4o  U • oZ 

tative  error  bars  for  the  data  are  shown.  The  values  of  T^  for  the  samples 
used  in  this  study  eure  listed  in  Table  I.  To  fit  the  "AlAs-like"  modes  we 
have  used  Eq.  (1)  and  (2)  with  F^^  =  2.8  cm  in  good  agreement  with  the 
300K  intrinsic  linewidth  of  the  end  point  material  AlAs.  The  solid  lines 
for  the  LOj^  modes  of  the  three  samples  in  Fig.  2  are  a  least  squares  fit  of 
Eqs.  (1)  and  (2)  to  the  experimental  data  using  only  the  "Spatial 

Correlation  Length"  of  this  mode,  as  an  adjustable  paurameter.  We  refer  to 
this  as  the  one  parameter  SCM  fit.  The  agreement  between  experiment  and  fit 
is  quite  good.  Values  of  (and  for  the  various  samples  used  in  this 

study  are  listed  in  Table  I. 

For  the  "InAs-like"  mode  of  sample  C  (T^=515°C)  the  data  also  can  be 
fit  by  the  SCM  using  F^^  =  2.8  cm  ^  (intrinsic  300K  linewidth  of  InAs)  and 
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varying  only  i.e.,  one  parameter  SCM  fit.  This  is  shown  by  the  solid 

line  in  Fig.  2  for  this  mode,  the  value  of  being  listed  in  Table  I.  For 

the  LO^  modes  of  the  other  samples  we  have  not  been  able  to  achieve  a  good 

fit  by  the  one  parameter  SCM  model.  For  these  features  we  have  had  to  vary 

not  only  but  also  rQ2-  This  two  parameter  SCM  fit  is  shown  by  the  dashed 

lines  for  LO  for  samples  A  and  F  in  Fig.  2.  The  values  of  F  .  and  1.  also 

2 

are  listed  in  Table  I  for  the  different  samples.  The  significance  of 

varying  in  addition  to  i  will  be  discussed  below. 

From  the  SCM  line-shape  fit  we  can  also  evaluate  the  full-width  at 

half-maximum  of  the  LO.  and  L0_  modes,  which  we  designate  at  P^,  and  P  , 

j.  2  tl  t2 

respectively.  These  values  also  are  listed  in  Table  I. 

For  the  LO^  mode  for  all  samples  except  C  the  peurameter  P^^  as  well  as 

2,2  had  to  be  varied  in  the  SCM  fit  (See  Fig.  2).  This  is  an  indication  that 

for  these  samples  a  new  scattering  mechanism  has  been  activated  by  the 

growth  conditions.  The  intrinsic  linewidth  P^  of  the  end  point  materials  at 

a  given  temperature  is  a  consequence  of  anharmonic  interactions  involving 

two  phonon  processes,  the  most  important  decay  channels  being  combinations 

23 

of  both  optical  and  acoustical  phonons  .  Reference  23  points  out  that  the 
intrinsic  phonon  lifetime  is  related  to  P^  through  the  energy-time  uncer¬ 
tainty  relationship.  Thus  a  change  in  P^  indicates  that  new  energy  relaxa¬ 
tion  channels  have  opened  up. 

In  addition  to-  the  RS  at  300°K  we  have  also  measured  in  PL  at  4°K. 

Figure  3  shows  the  dependence  of  the  PL  linewidth  on  T  .  All  the  samples 

13  ^ 

exhibited  the  same  peak  energy  .  The  values  of  the  PL  linewidth  also  are 

listed  in  Table  I  for  the  samples  of  this  study.  The  16  meV  linewidth  of 

sample  C  and  the  19  meV  linewidth  of  sample  D  are  significantly  narrower 

~  12  24 

than  previously  reported  values  (  25  meV)  for  this  material  ' 

Several  authors  have  discussed  the  origins  of  broadening  in  the  PL  of 

alloy  semiconductors.^^  Schubert  et  al.  have  investigated  the  effects  of 

thermal  broadening,  macroscopic  inhomogeneity  and  alloy  broadening  due  to 

random  cation  distribution  in  Ga,  Al  As^^.  References  15  and  16  have 

presented  statistical  models  based  on  this  random  distribution  to  account 

for  the  low  temperature  PL  linewidth,  the  two  calculations  arriving  at  the 

same  result.  In  Ref.  17  a  quantum  mechanical  formalism  was  used  to  evaluate 

the  exciton  volume,  thus  leading  to  more  realistic  values  of  the  exction 

15  16 

linewidth.  Singh  and  Bajaj  '  have  pointed  out  that  the  linewidth  is 
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strongly  dependent  on  the  short  range  disorder  present  in  the  alloy.  Raman 
scattering  provides  an  important  complement  to  PL  in  evaluating  alloy 
potential  fluctuations.  For  example,  the  PL  linewidth  may  also  be  influ¬ 
enced  by  impurities  while  RS  line-shapes  are  relatively  insensitive  to  this 

parameter.  The  LO  phonon  linewidth  from  the  depletion  layer  of  n-GaAs  for 

19  -3  25 

n=lxlO  cm  is  essentially  the  same  as  undoped  material. 

The  results  of  the  RS  and  PL  measurements  clearly  show  that  of  all  the 

material  studied  sample  C,  with  T^  =  515°C  exhibits  the  optimal  optical 

properties.  It  has  the  narrowest  PL  and  RS  linewidths,  has  the  largest 

and  in  addition,  has  the  smallest  (equal  to  that  of  the  end  point 

material)  so  that  the  one  parameter  SCM  can  be  used  to  fit  both  LO^  and  LO^ • 

Thus  we  conclude  that  under  this  growth  condition  there  is  a  decrease  in  the 

alloy  potential  fluctuations  of  this  material,  i.e.  a  tendency  towards 

ordering.  The  material  Al^  approximately  a  50-50  mixture 

of  AlAs  and  InAs  (AlInAs2)  and  hence  has  the  maximum  bond- length  and 

bond-angle  distortions  which  would  serve  to  drive  the  system  to  ordering 

3-5 

under  the  proper  thermodynamic  conditions 

It  has  been  suggested  recently  tha*-  the  asymmetry  in  the  LO  phonon 

line-shape  in  alloys  can  be  accounted  for  by  a  "Discrete-Continuum  Model", 

21  22 

i.e.  Fano-type  interaction  '  .  The  Fano  model  always  involves  two  parame¬ 

ters,  i.e.,  an  asymmetry  factor  (Q„)  and  broadening  parameter  (P  ).  In  Fig. 

r  F 

4  we  have  again  plotted  the  experimental  data  (dotted  lines)  for  samples  A, 
C  and  F.  The  solid  lines  for  both  LO^^  and  LO^  are  a  least-squares  fit  of 
the  data  to  Eq.  (2)  of  Ref.  22.  The  fit  using  the  Fano  model  is  quite  good 
and  is  comparable  to  the  SCM  fit  in  Fig.  2.  The  obtained  values  of  Q  and 

r 

r  for  these  samples  is  indicated  in  the  figure.  We  have  found  that  sample 

C  has  the  smallest  r„  for  both  modes.  The  main  result  of  this  comparison  is 

F 

to  demonstrate  that  as  far  as  fitting  procedures  are  concerned  the  two 
parameter  SCM  model  is  equivalent  to  the  two  parameter  Fano  model.  The 
authors  of  Ref.  22  have  compared  the  Fano  model  to  the  one  parameter  SCM 
without  explicity  pointing  out  that  the  former  intrinsically  involves  two 
variables  without  considering  a  two  parameter  SCM. 
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B.  Raman,  Transmission  Electron  Microscopy  and  Conductivity  Measurements  of 
Molecular  Beam  Deposited  Microcrystalline  Si  and  Ge. 

The  structure  of  molecular  beam  deposited  microcrystalline  silicon  and 
germanium  films  prepared  under  different  growth  conditions  has  been  analysed 
by  conventional  transmission  electron  microscope  (TEM)  and  Raman  scattering 
(RS)  and  compared  to  electrical  conductivity  experiments  .  The  TEM  mea¬ 
surements  yield  an  average  grain  size  ranging  from  20oS  to  1.5  nm.  On 
the  other  hand,  the  line  shape  of  the  RS  is  determined  by  a  mean  free  path 
related  to  the  average  separation  i.  between  defects  (or  impurities),  where 
I  15o8.  We  find  in  a  number  of  cases  that  the  electrical  conductivity  is 
determined  by  1  rather  than  L^.  This  experiment  demonstrates  the  signifance 
of  RS  as  a  structural  characterization  method  when  usedin  conjunction  with 
other  techniques  such  as  TEM. 

In  this  work.,  the  structure  of  molecular  beam  deposited  (MBD)  )ic-Si  and 

HC-Ge  films  grown  under  different  conditions  has  been  analyzed  by  TEM  and 

RS.  The  TEM  measurements  reveal  that  the  films  have  an  average  grain  size 

ranging  form  20o8  to  1.5^m.  In  contrast,  the  line  shape  of  the  Raman 

spectrum  (i.e.,  linewidth  and  asymmetric  broadening)  indicates  a  value  of  I 

much  smaller  than  the  grain  size  determined  by  TEM.  The  discrepancy  is 

attributed  to  generic  defects  or  impurities  introduced  during  growth.  These 

defects  (or  impurities)  destroy  the  translational  symmetry  of  the  crystal, 

reduce  the  mean  free  path  of  the  electron-phonon  system,  and  hence,  affect 

the  Raman  line  shape  in  a  manner  similar  to  that  observed  for  structural 

damage^^  or  finite  particle  size^®'^^.  We  find  that  in  our  films  il;£l5oR 

and  that,  in  a  number  of  cases,  it  dominates  the  electrical  conductivity. 

The  usual  q  =  0  selection  rule  of  first-order  Raman  scattering  is  a 

consequence  of  the  full  translational  symmetry  of  an  "ideal"  crystal^^.  The 

RS  of  an  "ideal"  diamond  or  zincblende  type  material  has  a  Lorentzian  line 

shape  with  and  intrinsic  linewidth  of  about  3  cm  ^  at  room  temperature. 

28  29  31 

Several  authors  '  '  have  reported  that  "finite  size  effects"  (such  as 

those  observed  in  uc-Si),  which  destroy  the  full  translational  symmetry  of 
the  material,  result  in  a  low-frequency  asymmetric  broadening  and  red  shift 

of  the  Raman  spectra;  however,  other  physical  phenomena  which  also  destroy 

27 

the  translational  symmetry  of  the  system,  such  a  structural  damage, 

19 

alloying,  etc.,  produce  a  similar  "q-vector  relaxation"  of  the  RS.  From 
the  details  RS  line  shape  a  parameter  can  be  determined  which  is  a  measure 
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of  the  average  xiistance  between  defects  (structural  damage)  or  finite 
particle  size  (nc-Si).  In  our  experiments  on  uc-Si  and  uc-Ge,  the  RS  line 
shape  is  determined  by  the  mean  free  path  of  the  electron-phonon  system 
induced  by  defects  (or  impurities)  rather  than  the  particle  size  effect.  On 
this  basis,  we  cannot  only  explain  why  the  determined  by  RS  is  much 
smaller  that  the  grain  size  evaluated  by  TEM,  but  can  arrive  at  a  quanti¬ 
tative  measure  of  the  mean  free  path.  In  general,  when  the  mean  free  path 
is  comparable  to  the  particle  size  L^,  the  fit  is  further  complicated  by  the 
Matthiessen ' s  rule,^^  where  most  of  our 

present  work,  thus,  the  Raman  spectra  are  mainly  influenced  by  1. 

It  has  been  shown  that  the  details  of  the  Raman  line  shape  (i.e.,  shift 

and  asymmetry)  associated  with  q-vector  relaxation  can  be  accounted  for  by  a 
...u _ _  -  _ _ 27,28  ^  1  _ _ 


model  which  assumes  a  Gaussian  spatial  correlation.' 


In  this  letter,  we 


extend  this  model  to  tcike  into  account  the  finite  instrumental  resolution  of 
the  Raman  system  in  order  to  make  a  more  accurate  comparison  with  experi¬ 
ment.  If  we  asssume  an  isotropic,  the  Raman  intensity  at  frequency  ui. 


I  ((ij),  can  be  written  as 
_  I 


where  q  is  expressed  in  units  of  2ii/a,  a  being  the  lattice  constant,  and 

r  =  3.0  cm  ^  is  the  intrinsic  Raman  linewidth  at  300K.  For  simplicity,  we 
o 

assume  a  triangular  instrumental  transfer  function  (the  integral  having  a 

33  - 1 

very  simple  analytical  form  )  with  a  full  width  at  half-maximum  r^(=lcm 

for  our  experimental  conditions).  For  the  dispersion  u)(q),  we  have  used  an 

analytical  form  for  a  one-dimensional  linear  chain  model  as  in  Eq.  (2). 

From  the  above  equation,  the  full  width  at  half-maximum  r(=  +r^)  and 

asymmetry  can  be  determined  as  a  function  of  I,  where  and  are 

the  low  and  high  frequency  half-widths,  respectively.  This  equation  differs 

from  Eq.  (1)  since  it  explicitly  includes  the  instrumental  transfer  function 

T  (^).  Thus  it  is  more  complete  than  Eq.  (1). 

Shown  by  the  solid  line  in  Figs.  5(a)  and  5(b)  are  F  and  3s  a 

function  of  I  for  the  case  of  Si  and  Ge,  respectively,  evaluated  from 
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Eq.  (3).  Also  shown  are  the  experimental  values  of  r  and  F  /F,  obtained 
r  a  b 

from  samples  prepared  under  various  depostion  conditions.  A  more  detailed 
description  of  sample  preparation  can  be  found  in  Ref.  34.  As  can  be  seen, 
the  agreement  between  experimental  and  calculated  F  and  values  for  Si 

and  Ge  is  quite  good.  Notice  that  for  all  the  measured  Si  and  Ge  samples, 
i£  i5o2. 

Figure  6  shows  the  TEM  micrograph  for  a  germanium  film  having  a  thick¬ 
ness  of  approximately  lOOoX.  The  film  was  deposited  at  a  substrate  tempera¬ 
ture  =  180°C  followed  by  thermal  annealing  at  550*^C  for  1  h.  From 
Fig.  6(b)  it  is  clear  that  that  average  grain  size  exceeds  lum.  In 
Fig.  6(a)  the  nearly  single  crystalline  electron  diffraction  patterns 
further  indicate  that  the  microcrystalline  grain  size  exceeds  the  electron 
beam  size;  however,  Raman  measurements  in  the  same  sample  show  that  the 
optical  phonon  branch  has  a  full  width  at  half-maximum  of  4.5cm  ^  and  an 
asymmetry  of  1.25.  According  to  Fig.  5(b),  this  corresponds  to  ll=15o8. 

The  comparison  of  these  two  experimental  results  shows  that  while  TEM 
reveals  the  dimensions  of  the  grain  structure,  RS  can  aibe  used  to  infer  the 
presence  of  local  structural  defects  through  a  parameter  dictated  by  1. 
Obviously,  the  detailed  microscopic  description  of  these  structural  defects 
is  not  embodied  in  our  approach. 

Figures  7  and  8  compare  room-temperature  electrical  conductivity  oand 
Raman  linewidth  F  measurements  for  wc-Si  and  uc-Ge  films,  respectively,  as  a 
function  of  isochronal  anneals.  The  data  points  were  obtained  from  samples 
prepared  at  various  T^.  For  silicon,  the  anneals  are  in  the  range  of 
600-800°C  while  for  germanium  the  range  varies  from  400  to  550°C;  in  the 
latter  material  the  lower  range  was  used  since  crystallization  in  germanium 
occurs  at  lower  temperatures  compared  to  Si  films. 

For  the  as-prepared  silicon  samples,  T^<360°C  leads  to  totally  amor¬ 
phous  structures,  at  380  T^  420°C  a  mixed  phase  of  amorphous  and  crystalline 
material  was  observed,  whereas  for  T^>420°C  a  complete  microcrystalline 
structure  is  obtained.  The  corresponding  temperatures  for  germanium  are 

T  <180°C,  180°C<T  <200°C,  and  T  >220°C,  respectively.  In  addition,  the 
s  s  s 

transition  from  amorphous  to  crystalline  phase  in  as-deposited  amorphous 

silicon  films  was  observed  at  annealing  temperatures  T  >600°C  for  one  hour 

depending  on  the  preparation  conditions.  For  amorphous  germanium  heat 

treatments  at  T  >450°C  lead  to  microcrystalline  samples. 

61 
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Indicated  in  Figs.  7  and  8  is  the  average  particle  size  (L  )  as  deter- 

o 

mined  by  TEM  analysis  after  the  annealing  800  and  550  C  respectively; 

however,  no  appreciable  change  in  was  observed  for  lower  annealing 

temperatures  as  long  as  no  amorphous  component  was  present  prior  to  the 

annealing.  Notice  that  the  values  for  L  range  from  20o8  to  1.5  urn  depend- 

34  ° 

ing  on  T  .  Tsu  et  al.  have  previously  used  a  quantitative  model  based  on 
s 

homogeneous  and  impurity  induced  heterogeneous  nucleation  to  explain  the 

wide  range  in  particle  size  observed  in  MBD-prepared  pc-Si  and  pc-Ge.  They 

found  that  the  presence  of  impurities  in  the  as-prepared  amorphous  films, 

particularly  at  lower  T  ,  results  in  smaller  particles  upon  crystallization. 

s 

The  top  parts  of  Figs.  7  and  8  show  o  for  Si  and  Ge  films  respectively  as  a 

function  of  the  annealing  temperature.  For  Si,  after  annealing  at  800°C,  a 

total  increase  in  of  approximately  two  orders  of  magnitude  was  observed 

for  samples  having  L^''l.5Mm.  The  total  increase  in  conductivity  gradually 

decteases  as  the  average  particle  size  decreases.  In  samples  prepared 

either  at  T  >600°C  or  at  T  <200°C  no  significemt  change  in  o  was  observed, 
s  s 

The  increased  in  for  intermediate  T^  is  attributed  to  the  annealing  of  some 
of  the  structural  defects  introduced  during  the  crystallization  process 
occurring  at  lower  temperatures.  Raman  measurements  [Fig.  7(b)]  on  the  same 
samples  after  each  annealing  step  correlate  well  with  conductivity  measure¬ 
ments.  For  intermediate  T  ,  the  Raman  linewidth  T  decreases  as  T  is 

s  a 

increased.  According  to  Fig.  5(a),  the  decrease  in  F  corresponds  to  an 
increase  in  the  length  £.  The  similar  trend  of  a  and  I  indicates  a  reduc¬ 
tion  in  the  defect  density  after  annealing. 

Previously  we  found  that  the  void  structure  of  MBD-prepared  silicon 
films  deposited  at  Tg>350°C  is  much  reduced,  therefore  limiting  the 
indiffusion  of  impurities.  Thus,  we  believe  that  the  increase  in  both  and 
is  related  to  defects  rather  than  impurities.  At  high  deposition  tempera¬ 
ture,  for  instance,  T  =  600°C  in  Fig.  7,  the  structural  defects  are  an- 

s 

nealed  during  the  growth  process  and  therefore  no  change  in  either  conduc¬ 
tivity  or  Raman  linewidth  is  observed  for  post-deposition  heat  treatments. 
The  sitution  is  quite  different  for  samples  having  small  particles.  At  T^  = 
200°C,  the  particle  size  is  in  the  range  of  20o8  as  detemined  by  TEM,  Heat 
treatments  do  not  change  significantly  either  the  electrical  conductivity  or 
the  Raman  linewidth.  What  this  suggests  is  that,  unlike  samples  having 
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larger  particle _si2e,  the  dominant  scattering  mechanism  is  due  to  grain 
boundaries . 

For  Ge,  in  the  range  of  annealing  temperatures  used  as  shown  in 
Figs.  8(a)  and- 8(b),  no  apprecicible  change  in  either  a  or  r  was  observed. 
This  may  be  due  to  the  fact  that  there  are  fewer  generic  defects  present  in 
Ge  as  compared  to  Si  although  another  possibility  is  that  the  structural 
defects  in  germanium  are  more  stcible  than  in  silicon  and  perhaps  annealing 
temperatures  in  excess  of  550°C  are  needed  to  observe  the  thermal  annealing 
of  the  defects. 

C.  Symmetry  Forbidden  LO-phonon  Raman  Scattering  from  Electric  Fields  in 
the  Space  Charge  Region  in  Heavily  Doped  <100>  n-GaAs. 

We  have  investigated  symmetry  forbidden  longitudinal  optic  Raman 
scattering  for  various  polarization  conf iguratioins  at  4579A  from  the  <100> 
surface  of  heavily  doped  n-GaAs  caused  by  the  strong  surface-fields  and  high 
impurity  levels^^.  We  have  evaluated  the  magnitude  and  phase  of  the  coeffi¬ 
cient  of  the  electric-field  induced  term  as  well  as  the  magnitude  of  the 
coefficient  of  the  impurity- induced  factor.  These  parameters  can  be  used  to 
gain  information  edxjut  the  surface  electric  fields  in  <100>  n-GaAs. 

We  have  studied  the  symmetry  forbidden  LO  phonon  RS  from  the  <100> 
surface  of  heavily  doped  n-GaAs: Si  due  to  the  high  impurity  levels  and 
strong  surface  electric  fields.  The  samples  were  grown  by  molecular  beam 

epitaxy  and  had  free  carrier  concentration  (electrons)  in  the  range 
18  “3  19  *“3 

1  X  10  cm  <n<l  X  10  cm  .  The  free  carrier  concentration  as  well  as 

space  charge  region  (SCR)  electric-fields  have  been  well  characterized  in 
36  37 

previous  studies  '  .  The  measurements  were  made  at  room  temperature  in 

the  backscattering  geometry  using  the  457 9A  line  of  an  Ar-ion  laser.  We 
have  studied  RS  for  all  four  possible  polarization  configurations,  i.e., 
I^[x(y,z)x],  l2[x(y,y)x),  l2lx(y' ,y' )x) ,  and  I^lx(z' ,z' )x] ,  where  x  =  <100>, 
y  =  <010>,  z  =  <001>,  y'  =  <011>  and  z'  =  <011>.  From  these  measurements  we 
have  evaluated  the  coefficients  of  the  electric  field-  and  impurity- induced 
terms  at  4579  A. 

Shown  in  Figure  9  is  the  RS  from  sample  782-8A  in  the  normally  allowed 
configuration  Ij^[x(y ,2)x] .  The  peaks  designated  and  L_  are  the  coupled 

plasmon  -LO  phonon  modes  which  originate  from  the  bulk  of  the  material  where 
the  free  carriers  exist.  The  structure  marked  LO  (292  cm  ^)  is  the 


unscreened  tiO  phonon  mode  of  GaAs  which  originates  only  from  the  SCR.  In 

previous  studies  the  L  ,  L_  modes  have  been  used  to  characterize  the  free 

36  37  ~ 

carrier  concentration  '  (and  hence  impurity  concentration).  Since  GaAs 

18  -3 

is  degenerate  for  n>lxlO  cm  we  take  the  impurity  concentration,  n‘ ,  equal 
to  the  free  carrier  concentration,  n.  The  amplitude  of  LO(292  cm”  ),  in 
relation  to  that  of  an  undoped  sample,  has  been  utilized  to  evaluate  the 
width  of  the  SCR  and  hence  the  built-in  surface  electrilc  field^^.  In  this 
study  we  will  concentrate  on  the  polarization  dependence  of  the  unscreened 
LO  (292  cm  mode. 

In  Fig.  10  we  display,  by  the  points,  the  experimentally  determined 

ratios  1^/1^,  ^  function  of  impurity  concentration  n^  for 

a  laser  excitation  wavelength  of  4579A.  Note  that  as  n.-KD,  the  ratio 

12/1^^'*^  while  both  1^/1^  and  We  will  show  below  that  difference 

between  13/1^  and  with  increasing  n^  is  due  to  an  interference  effect 

similar  to  that  reported  in  Refs.  38  and  39. 

The  effects  of  electric  fields  and  linear-q  terms  on  the  Raman  tensor 

40  41 

have  been  discussed  in  several  articles  '  .  We  will  demonstrate  below 


that  linear-q  terms  are  not  significant  in  our  experiment  and  hence  we  will 
consider  only  the  electric  field-induced  effects.  Impurity- induced  terms 


will  be  discussed  in  the  next  paragraph.  For  the  case  of  LO-phonon  scatter¬ 
ing  from  <100>,  with  the  electric  field  ^(x)  also  along  this 
crystallographic  direction,  the  Raman  tensor  can  be  written  as"^°'^^: 


In  Eq.  (4)  a^  (a2)  and  (p^  (^2)  are  the  magnitude  and  phase,  respectively,  of 
the  electric  field-induced  term  whild  d  and  are  related  to  the  normally 
allowed  atomic-deformation  (AD)  scattering  process.  For  our  scattering 
geometry,  i.e.,  from  the  <100>  surface,  it  can  be  shown  from  Eq.  (4)  that 
only  a^,  <p2»  and  can  be  observed. 
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as 


The  Raman  tensor  for  impurity- induced  scattering  can  be  written 
38,39, 


/ 

\ 


1 

0 

0 


0  0 


1  0 


o  1/2  ,  .  , 

Cn^'  exp{i<p^) 


0  1 


(5) 


where  C  and  (p^  are  the  magnitude  and  phase,  respectively,  of  the  coefficient 

for  this  process.  It  has  been  demonstrated  by  Menendez  and  Cardona  that 

38  38 

Eqs.  (4)  and  (5)  represent  independent  scattering  processes  '  .  There¬ 

fore,  in  order  to  calculate  the  total  scattering  intensity  the  two  processes 
must  be  added  after  squaring,  not  before  squaring. 

From  Eqs.  (4)  and  (5)  we  can  write  expressions  for  the  ratios  of  the 
intensities  of  the  various  polarization  configurations; 


(6a) 


(6b) 


(6c) 


o 
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where  (p  =  (^>  -  (p  ) ,  L  is  the  width  of  the  depletion  region  and  a(X)  is  the 

^  -J  O 

penetration  depth  of  the  light  at  wavelength  X,  respectively.  In 

Eqs.  (6)  we  have  neglected  the  small  difference  in  wavelength  between  the 

incident  and  scattered  light.  Expressions  for  L  and  S(x)  as  a  function  of 

O 

n.  are  given  in  Ref.  36.  The  quantity  a(X)  is  obtained  from  the  work  of 
^  42 

Aspnes  and  Studna  .  Note  that  Eqs.  (6b)  and  (6c)  exhibit  an  interference 
between  the  electric-field- induced  and  AD  terms  similar  to  that  discussed  in 
Refs.  38  and  39. 

The  solid  lines  in  Fig.  10  are  at  least  squares  fit  of  Eqs.  (6)  to  the 
experimental  points.  From  such  and  analysis  we  can  evaluate 


a^/d  =  (5.0  ±  1.0)  X  lO"^  cm/V 
<p  =  60  ±  10° 

P  2  -20  3 

C  /d  =  (1.4  ±  0,2)  X  10  cm 

43 

Grimsditch  et  al  have  measured  d  and  also  have  derived  a  theoretical 
expression  from  which  both  d  and  cp^  can  be  evaluated  as  a  function  of  X. 
For  4679A  both  experiment  and  theory  yield  d  =  200A^  while  =  -2°  is 
obtained  from  theory.  Thus  we  have  that: 


a^  =  (9.6  ±  0.4)  X  10  ^  cm-A^/V 

<p  =  62  ±  10° 

=  (5.6  ±  0.8)  X  10'^°  cm^  -a"^ 

It  has  been  shown  that  forbidden  LO-scattering  induced  by  an  electric 

44 

field  can  be  written  as: 


a^ei‘p2-g3/g  ^3  Xij(u))l 


(7) 
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where  "Ka)  is' the  photon  energy  and  is  the  Linear  susceptibility.  In 

our  case  we  are  at  an  energy  near  the  E.  optical  transition  of  GaAs. 

45  ^ 

Trommer  et  al  have  used  a  two-dimens ioonal  critical  point  model  for  the 
transitions  so--  that 


3^/a 


[w 


(oj-oj^-in) 


-3 


(8) 


where is  the  energy  of  the  transition  in  GaAs  (-2.89  eV)  and  n  is  a 
broadening  parameter  taken  to  be  70  meV.  Thus,  from  Eqs.  (7)  and  (8)  we  can 
calculate  at  4579A  (=2.71  eV)  to  be 


<p2  (4579  A)  =  64° 


This  value  of  cp^  is  in  good  agreement  with  our  experimental  value  of 
62  ±  10°. 

Trommer  et  al^^  have  evaluated  a  /d  =  1.2  x  10  cm/V  at  4416A 

2  ^ 

(2.81  eV).  From  Ref.  43  d  =  280  A  and  hence  we  can  deduce  that 

-4  2 

(4416A)  =  3.4  x  10  cm-A  /V.  From  Eqs.  (7)  and  (8)  we  can  determine  that 

the  ratio  a_  (4416A)/a-(4579A)=6,  in  reasonable  agreement  with  the  experi- 

^  ^  -5 

mental  ratios  of  (3.4  x  10  /9.6  x  10  )=3.5. 

Comparison  of  the  quantity  C  with  other  experiments  is  somewhat  more 

2  46  -12  -3  4 

difficult.  From  Ref.  39  C  can  be  evaluated  to  be  1.4  x  10  cm  -A  at 

1.85  eV.  This  value  is  considerably  larger  than  that  obtained  from  our 

experiment,  i.e.  5.6  x  10  ^^cm  ^-A^.  This  is  not  too  surprising  since  C^(w) 

has  maxima  under  resonance  conditions  and  the  value  of  Ref.  8  was  deduced 

from  resonance  with  the  E  +  L  optical  transition. 

o  o 

We  have  neglected  linear-q  terms  based  on  the  results  of  a  study  of  the 

27 

Raman  spectra  of  <100>  GaAs  implanted  with  various  fluences  of  As-ions 
It  was  found  that  the  highest  fluences  produced  and  undamaged  region  of 
dimensions  45-50A  but  that  no  violation  of  the  poleurization  selection  rules 
for  this  fluence  was  observed.  Therefore,  we  conclude  that  for  dimensions 
of  the  order  finite  size  effects,  i.e.,  linear  -q  terms  are  not  significant. 
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For  n.  =  l  X  -10  cm  ,  the  highest  impurity  level  used  in  this  study,  the 
width  of  the  SCR  is  about  50A  .  Therefore,  we  have  neglected  linear  -q 
terms  in  our  analysis. 

D.  Comprehensive  Investigation  of  Polish-Induced  Surface  Strain  in 
<100>  and  <111>  GaAs  and  InP 

We  have  performed  a  systematic  and  extensive  investigation  of 

polish-induced  surface  strain  in  <100>  and  <111>  GaAs  and  InP  using  Raman 

47-49 

scattering  from  the  Longitudinal  optical  (LO)  phonon  modes  .  By  using 

various  lines  of  an  Ar-ion  laser  it  was  possible  to  accomplish 

47-49 

non-destructive  depth  profiling  .  To  account  for  the  observed  lineshape 
cheuiges  we  have  used  a  mode  which  is  based  on  the  convolution  of  the  pene¬ 
tration  depth  of  the  light  and  skin  depth  of  the  polish- induced  surface 
strain.  From  such  an  analysis,  we  have  obtained  the  polish  induced  surface 
strain,  skin  depth  of  the  strain  and  inhomogeneous  broadening.  For  the 
<100>  surface,  the  strain  is  about  2-3%  in  both  materials  and  the  skin  depth 
(  100-500A)  is  relatively  independent  of  particle  size.  In  contrast,  for  the 
<111>  surface,  the  average  surface  strain  is  only  about  0.6%  for  GaAs  and 
1.2%  for  InP  and  the  skin  depth  is  of  the  order  the  particle  size.  The 
dependence  of  the  strain  skin  depth  on  polish  time  also  has  been  studied.  A 
qualitative  argument  based  on  polish- induced  bond  breaking  is  proposed  to 
explain  why  surface  strain  for  <111>  is  considerably  less  than  for  <100>  and 
why  the  strain  is  compressive  for  both  surfaces.  Using  a  one  dimensional 
diffusion  model,  we  can  successfully  explain  the  depth-dependence  of  the 
polish- induced  strain  and  the  polish  time  dependence  of  the  damage 
skin-depth.  This  analysis  yields  a  diffusion  coefficient  for  the 
polish- induced  strain  for  the  two  surfaces.  Our  diffusion  model  is  consis¬ 
tent  with  the  conventional  model  of  chemo-mechanical  polishing  of  compound 
semiconductor . 

Polish- induced  surface  damage^^  in  diamond  and  zincblende-type  semicon¬ 
ductors  has  been  investigated  by  several  methods  including  optical  tech- 

51  52-54  55 

niques  such  as  photo luminescence  ,  reflectance  ,  absorption  ,  photon 

backscattering^^,  double  crystal  X-ray  diffractometry  ,  transmission 

57  57'”59 

electron  microscopy (TEM)  ,and  Raman  spectroscopy  (RS)  .  The  technique 

of  RS  is  a  extremely  useful  for  a  number  of  reasons.  The  influence  on  the 

40,60  , 

RS  spectrum  of  damage- induced  perturbations  such  as  strain  and 
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disorder  (which  produce  different  effects)  are  well  understood.  Nonde¬ 
structive  depth  profiling  can  be  done  using  various  laser  lines.  Since  the 
Longitudinal  Optical (LO)  phonon  in  zincblende-type  materials  is  a  singlet, 
certain  complications  which  occur  in  electronic  transition  spectroscopies 
are  avoided  such  as  strain- induced  interband  and  intraband  splitting^^. 

1.  Measurement  Conditions 

The  method  of  preparation  of  the  samples  is  described  in  Ref.  47.  All 
measurements  were  made  at  room  temperature,  using  an  Ar-ion  laser  as  an 
excitation  source  with  output  wavelengths  at  5145,  5017,  4880,  4765,  and 
4579  A.  In  order  to  avoid  heating  effects  a  line  focus  (300|im*3mm)  was  used 
in  which  the  power  was  300mW  (except  for  4579  and  5017  A,  in  which  case  it 
was  130mW).  For  the  <100>  surface,  the  allowed  cross  polarization  configu¬ 
ration  [x(y,z)x]  was  employed  to  reduce  the  Rayleigh  scattering  from  the 
rough  surface,  where  x,y,z  indicate  <100>,  <010>,  and  <001>,  respectively. 
The  cross  configuration  for  the  <111>  surface  is  forbidden  so  it  was  not 
analyzed.  We  found  no  significant  difference  in  the  line  shape  for  differ¬ 
ent  polarizations. 

2.  <100>  Surface 

Shown  in  the  Fig.  11(a)  is  the  RS  using  the  5145  A  line( labeled  A)  in 
the  region  of  LO  phonon  for  undamaged  <100>  InP  and  material  polished  with 
0.05  pm  grit  for  the  three  different  times  (5,  10,  and  60  minutes).  The 
undamaged  RS  is  fairly  shau^),  with  a  symmetry  peak  center  at  345.6  cm 
Spectrum  B,  for  which  the  polishing  time  is  5  minutes,  has  a  relatively 
sharp  feature  at  about  348  cm  ^  and  a  broad  shoulder  at  eibout  355  cm  ^ .  As 
the  polish  time  increases,  the  amplitude  of  the  shoulder  grows  relative  to 
the  low  frequency  feature,  but  the  position  of  the  shoulder  does  not  shift. 

The  spectra  in  Fig.  11(b)  are  the  results  for  <100>  GaAs.  Spectrum  A  is 
for  undamaged  material  taken  with  a  5145A  line,  and  spectra  B-D  are  for  the 
case  of  0.05  urn  grit  polish  for  three  different  times  (5,  10,  60  minutes, 
respectively).  The  undamaged  spectrum  has  a  sharp  peak  at  292.1  cm  .  We  see 
that  with  increasing  polish  time,  the  amplitude  of  the  shoulder  grows 
relative  to  the  low-frequency  feature,  and  the  position  of  the  shoulder 
shifts  to  higher  frequencies.  We  will  demonstrate  later  in  the  paper  that 
the  two  features  in  the  spectra  of  damaged  InP  and  GaAs  do  not  arise  from 
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two  separate  peaks  but  are  the  consequence  of  the  convolution  of  the  skin 
depth  of  the  damage- induced  strain  and  the  optical  penetration  depth  for  the 
various  wavelengths. 

Plotted  in  Fig.  12  are  the  results  for  <100>  InP  and  GaAs  for  fixed 
time  and  different  wavelengths.  The  InP  was  polished  with  0.05  4m  grit  for 
5  minutes  as  shown  in  Fig.  12(a),  while  GaAs  was  polished  with  0.3  4m  grit 
for  60  minutes  as  shown  in  Fig  12(b).  The  spectra  were  taken  with  5145, 
5017,  4880,  4765  and  4579  A  excitation  lines,  respectively,  for  both  materi¬ 
als.  The  amplitude  of  RS  of  the  damaged  material  was  consideraibly  less  than 
the  undamaged  sample  although  the  integrated  intensities  remained  approxi¬ 
mately  constant.  We  can  see  that  as  the  penetration  depth  of  the  light 
decreases,  the  amplitude  of  the  high  frequency  shoulder  grows  relative  to 
the  low  frequency  feature.  For  GaAs  the  high-frequency  shoulder  in  spectrum 
B  has  become  the  dominant  feature  at  the  shortest  wavelength. 

3.  <111>  Surface 

Plotted  in  Fig  13(a)  are  the  RS  in  the  region  of  the  LO  phonon  for 
undamaged  <111>  InP  using  the  5145-A  line  (spectrum  A)  and  polished  with 
0.054m  grit  for  10  minutes  using  5145,  4880,  and  4579  A  excitation  lines 
(labeled  B,  C,  and  D  , respectively ) .  We  have  not  plotted  the  spectra  for  all 
wavelengths  since  there  are  no  significant  differences  between  them.  The 
<111>  and  <100>  [see  Fig.  12(a)]  undamaged  RS  are  identical  in  terms  of  peak 
position  and  linewidth.  However,  for  damaged  conditions,  the  RS  from  the 
<111>  surface  are  quite  different  from  those  from  the  <100>  surface.  The 
relatively  sharp  feature  at  about  348  cm  which  appeared  in  the  RS  from 
the  <100>  surface  [  see  Fig  12(a)  ],  is  not  seen.  The  broadening  is  fairly 
symmetrical,  and  spectra  peak  at  about  346.5-347.5  cm  As  the  penetration 
depth  of  the  light  decreases,  the  linewidth  increases  and  the  peak  position 
shifts  to  the  blue. 

The  RS  spectra  in  the  region  of  the  LO  phonon  for  undamaged  and  damaged 
<111>  GaAs,  polished  with  0.3  4m  grit  for  10  minutes,  are  shown  in  Fig 
3(b). As  for  case  of  InP,  the  undamaged  spectra  for  <111>  and  <100>  surfaces 
are  identical.  Also  for  GaAs  the  damaged  spectra  are  symmetrically  broad¬ 
ened;  with  decreasing  wavelength  the  RS  features  broaden  and  shift  to  the 
blue,  although  not  as  much  as  for  InP  [see  Fig  13  (a)]. 
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The  measured  RS  shift  and  line-width  of  the  <111>  InP  and  GaAs  for 

various  polishing  conditions  versus  different  penetration  depths  of 

the  different  Ar-ion  laser  wavelengths  (X)  are  plotted  in  Fig.  14(a)  and 

14(b),  respectively.  The  RS  shift  and  linewidth  for  <111>  InP  and  GaAs 

damaged  with  0.3^m  grit  for  60  minutes  do  not  change  with  different  d  . 

o 

However,  for  the  O.OSijm  grit  size  the  RS  shift  and  linewidth  do  increase 
considereibly  with  decreasing  light  penetration  depth.  For  the  undamaged 
surfaces,  the  RS  shift  and  linewidth  are  constant  with  d^,  demonstrating 
that  a  strain-free  surface  can  be  achieved  by  proper  mechanical  and  chemi¬ 
cal  polishing. 


4.  Analysis 
a.  Convolution  Model 

It  is  well  known  that  polishing  induces  strain  at  the  surface  of  a 
material  thus  producing  a  shift  in  LO  frequency  from  its  unstrained  val- 
^^40,47  49,60,^  damaged- induced  strain  (and  hence  RS  shift)  diminishes 
as  a  function  of  depth  from  the  surface.  The  observed  RS  signal  is  the 
convolution  of  the  depth  dependence  of  the  LO  frequency  and  the  optical  skin 
depth  (d^)  of  the  light.  Assuming  an  exponential  depth  dependence  for  the 
Rameui  frequency 


a)(x)=u)  +  Aw  exp(-x/d  ) 

O  0*^  8 

(9) 

w(x)=  Aw^exp( -x/d^) 

(10) 

where  x  is  the  dist2Uice  from  the  surface,  is  the  unperturbed  frequency, 

Ao)  is  the  shift  due  to  surface  strain  and  d  is  the  skin  depth  of  the 
o  s 

strain.  The  Raman  intensity  I(u))  is  given  by 


I((jj)a 


dxL[uj-u(x)  ]exp(-2x/d^) 


(11) 


!> 

where  L((ij)  is  a  line  shape  factor  which  is  a  convolution  of  the  lorentzian 
line  shape  of  the  RS  and  the  instmimental  transfer  function.  We  have  ne¬ 
glected  the  small  difference  in  the  optical  skin  depth  for  the  incident  and 
scarttered  radiation  in  the  above  discussion.  After  changing  variables  we 
write: 
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(12) 


<i)(Y)=u>  Y 

e  O  O 


where  y=exp(-x/d^) .  Equation  (11)  then  becomes; 


.00 


I(u))  a 


dYL[oj-(jj(Y)  1y 


P-1 


(13) 


where  p=2d^/d^.  For  simplicity  we  have  assume  a  triangular  instrumental 
transfer  function  T(5)  with  full  width  at  half  maximum  of  F^.  The  quantity 
L(u))  can  be  written  as: 


L{a))  = 


J  Ci^~  £)V(%)^ 


(14) 


where  F  is  the  linewidth  of  the  RS.  Equation  (14)  can  be  solved  to  yield  the 
analytical  form 


L(a))  =  (F/2rtrj)(R(a)+Fj)>R(u)-rj)-2R(u))]  (15) 

where 

R(u))  =  (2aj/r)arctan(2cj/r)-(l/2)ln[(2w/r)^+l]  (16) 


After  integration  by  parts  we  obtain: 

I(oj)a 


(17) 


where 


I^=R(4>+rj)+R(<t>-Fj)-2R(<l>)  (18a) 

<t=  uj-(aj  +A<ij  )  (18b) 

'  o  o' 

l2=A{arctan[ 2 ( H+F^ ) /F ] +arctanl 2 ( O-F ^ ) /F ] 

-2arctan(  2<>/F) }  (18c) 


with 


A=2Aoj^/r(li+l) 


{ 19a) 
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(19b) 


l3=B  dYY‘^''^{l/H  +  (6+o)^]+i/[i+(6-o)^]-2/(l+6^)} 


B=4(Aw^)^/r^((3+l) 

(20a) 

5=2[u)-(w^+A(jj^Y)  ]/r 

(20b) 

a=2rj/r 

(20c) 

It  is  also  possible  to  consider  any  inhomogeneous  strain  in  the  plane 
perpendicular  to  the  polishing  direction  (for  a  given  value  of  depth).  By 
increasing  the  broadening  parameter  r  in  Eg. (14)  from  its  intrinsic  value 
(2.4  cm  ^  for  GaAs  and  1.2  cm  ^  for  InP)  and  assuming  a  depth  dependence 
similar  to  that  of  the  frequency  shift  [see  Eg. (9)],  we  can  write; 

r=r  -*-Ar^exp(-x/d^)  (2i) 

where  is  considered  to  be  the  inhomogeneous  broadening. 

Two  dimensional  strain  in  the  plane  perpendicular  to  the  polishing 

direction  is  assumed  in  order  to  correlate  the  RS  shift  Aw  with  the  surface 

o 

strain.  The  quantity  for  <100>  or  <111>  surface  is  related  to  the 

surface  strain,  S,  by: 


where 


^100~^12'^^11 


(22b) 


and 


(  Au)^/w^)^^^={[(p+2q)(2-A^^^)6w^^)-[2r(UA^^^)/3u)^^]}S=K^^^S  (22c) 

where 


(22d) 
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In  the  above  equations,  C, ,  ,  C,_,  and  C,.  are  elastic  stiffness  constants 
^  11  12  44 

and  p,  q,  and  r  are  coefficients  which  describe  the  strain-induced  shift  of 

40 

the  optical  phonons.  For  GaAs  ^111~”*^’^  ’  parameters  p,q, 

and  r  have  not  been  determined  for  InP;  but  since  these  parameters  do  not 

40 

vary  much  for  the  III-V  materials  ,  we  use  GaAs  values  for  InP  and  hence 
K^qq=-1.3  and  this  semiconductor, 

b.  Fitting  of  <100>  Surface  Results 

Shown  in  Fig  15(a)  is  a  least-square  fit  of  Eq(9)-(20)  to  the  spectra 
at  5145  A  and  4880  A  of  <100>  InP  damage  with  0.3  urn  grit  for  lOmin.  The 

adjustable  parameters  were  Au)  {  RS  shift  due  to  surface  strain  )  and  d  .  We 

-1  °  -1  -1  ^ 
have  taken  u)^=345.6  cm  and  r^=1.2  cm  ;  r^=2.4  cm  from  slit  settings. 

The  d^  value  is  from  Ref.  42.  There  is  excellent  agreement  between  theory 

and  experiment  at  both  5145  and  4880A  wavelengths.  We  find  that  A  w  =12.3 

-1  -1  ° 
cm  and  d  =386A  for  5145A  and  A  w  =12.3  cm  and  d  =380A  for  4880A.  Listed 
s  os 

in  Table  II  are  summaries  of  the  results  for  InP  (averaged  over  the  differ¬ 
ent  wavelengths)  for  the  various  conditions. 

In  Fig.  15(b)  are  plotted  experiment  and  fit  for  <100>  GaAs  at  5145A 

and  4880A  damaged  with  0.3  pm  grit  for  10  mins  using  w  =292.1  cm  F  =2.4 
-1-1  o  o 

cm  ,  and  r_  =2.4  cm  .  It  is  found  that  A  w  =10.7  cm-1  and  d  =235A  for 
I  .  o  s 

5145A  and  Aw  =  10.9  cm  and  d  =  240  A  for  4880  A.  The  data  for  <100> 
o  s 

GaAs  for  different  polishing  conditions  are  listed  in  Table  II. 

An  inhomogeneous  strain  for  <100>  for  both  GaAs  and  InP  was  also  intro¬ 
duced.  We  find  that  increasing  AP^  to  3  cm  ^  had  little  effect  on  the 
quality  of  the  fits  in  Figs.  15(a)  and  15(b);  further  increases,  however, 
result  in  a  significant  deviation  in  the  correlation  between  fit  and  experi¬ 
ment.  Therefore,  we  conclude  that  the  inhomogeneous  strain  is  less  than  3 
cm  ^  for  the  <100>  situation.  It  is  significantly  less  than  the  homogeneous 

shift.  Aw  ,  which  is  8-12  cm  ^  (see  Table  II). 
o 

For  the  case  of  0.05  pm  grit  we  have  performed  a  systematic  study  of 

the  effect  of  the  polishing  time,  t,  on  A  as  shown  in  Table  II.  Note 

that  while  and  A  w  and  Aw  are  almost  independent  of  polishing  time 

o  o 

(greater  than  2.5  minutes)  the  quantity  d^  saturates  with  increasing  t. 
Plotted  in  Fig. 16  by  the  open  circles  are  the  experimental  values  of  d^  as  a 
function  of  t  for  <100>  polishing. 
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c.  Fitting  of  <111>  Surface  Results 

The  interpretation  of  the  experimental  results  for  the  <111>  case  are 

somewhat  different  from  that  for  the  <100>  situation.  We  were  not  able  to 

obtain  a  good  fit  to  the  experimental  data  if  we  used  only  Eq.(17),  i.e.  no 

homogeneous  broadening.  Therefore  it  was  necessary  to  use  Eq(9).  For  this 

case  inhomogeneous  broadening  is  an  important  effect,  as  found  in  the 
58 

previous  RS  study  .  In  contrast  to  the  <100>  case  it  is  difficult  to  get 
meaningful  data  from  the  spectrum  at  just  a  single  wavelength  since  varia¬ 
tions  in  d^  and  r  produce  similar  effects  on  the  theoretical  fit.  Therefore 

we  have  obtained  values  of  Aoj  ,  AF  and  d  for  the  different  polishing 

o  o  s 

conditions  by  performing  a  self-consistent  least  squares  fit  of  Eqs(l)-(11) 
to  all  five  spectra  at  different  wavelengths,  i.e.,  5145,  4965,  4880,  4765, 
and  4579  A.  The  resultant  numbers  are  listed  in  Table  III. 

Comparing  Tables  II  and  III,  we  find  that  for  the  <111>  surface,  the 
compressive  surface  strain  is  less  than  that  in  the  <100>  case,  the  dominant 
strain  is  inhomogeneous,  and  the  depth  of  damage -induced  strain  is  much 
deeper  than  that  of  the  <100>  surface  for  corresponding  polishing  condi¬ 
tions. 

For  this  surface  also  we  have  performed  a  systematic  study  of  the 

effects  of  polishing  time  on  Au,  AF  and  d  for  O.OSum  grit.  As  for  the 

o  s 

<100>  case  Aoi  ,  AF  are  independent  of  t  (  for  t>  2.5  mintues)  while  d  is 
o  o  s 

a  function  of  t  (see  Table  III).  The  dependence  of  d^  on  t  for  the  case  of 
the  <111>  surface  and  0.05pm  grit  is  illustrated  in  Fig.  17  by  the  open 
circles. 

5.  Discussion 

Our  results  show  that  the  surface  strain  is  compressive  and  is  about 
2.2%  for  <100>  GaAs,  2.6%  for  <100>  InP,  0.6%  for  <111>  GaAs  and  1.2%  for 
<111>  InP.  This  was  found  from  the  sign  and  magnitude  of  Aoo^  in  Tables  II 
and  III  and  Eqs.(ll)  and  (12),  The  inhomogeneous  strain  (  AF^/2)  is  less 
than  0.3%  for  the  <100>  surface  while  it  is  about  1.4%  for  the  <11 1>  sur¬ 
face.  The  strain  skin  depth  is  dependent  on  polish  time  and  it  approaches  a 
limiting  value  for  a  long  time  polish  with  a  given  grit  size. 

We  propose  an  ideal  spring  model  to  qualitatively  account  for  the 
compressive  nature  of  Aoi^  and  the  difference  in  Acj^  between  the  two 
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surfaces.  In  addition,  in  order  to  explain  the  saturation  and  time  depen¬ 
dence  of  d^,  we  propose  a  model  based  on  one  dimensional  diffusion  of  the 
strain. 

We  also  considered  the  effect  on  the  LO  phonon  line  shape  of  other 
perturbations  such  as  disorder,  temperature  gradient  effects  and  free 
carrier  interactions (coupled  plasmon-LO  modes),  a.  Ideal  Spring  Model 

Shown  in  Figs.  18  and  19  are  the  zincblende  crystal  structure  in  the 
<100>  and  <111>  directions,  respectively.  For  the  former  case  we  display  six 
layers,  from  n  to  n+5,  while  for  the  latter  situation  eight  layers  are 
shown,  from  n  to  n+7.  The  plane  designated  n  in  both  figures  is  not  neces¬ 
sarily  the  top  plane,  there  being  n-1  layers  above  it.  The  bond  distance  is 
(r3/4)  a  ,  where  a  is  the  lattice  constant.  In  discussion  below  we  assume 
that  the  bonds  between  atoms  are  ideal  springs  and  that  the  effect  of  the 
polishing  is  to  breeik  the  bonds  between  two  given  sets  of  layers  shown  in 
these  figures. 

For  the  <1Q0>  case  ail  layer  are  equivalent,  the  atoms  between  each 
successive  plane  being  connected  by  two  bonds.  We  assume  that  the  polishing 
remove  all  the  layers  above  n+1,  i.e.  the  bonds  between  n  and  n+1  are 
broken.  The  layer  n+1  will  then  be  pulled  toward  the  n+2  plane.  If  this 
displacement  is  complete,  i.  e.,  the  n+1  layer  is  pulled  entirely  into  the 
n+2  plane,  the  bond  distance  between  n+1  layer  atoms  and  n+2  layer  atoms 
become  {r2/4)  a  .  This  corresponds  to  a  compressive  strain  of  18%. 

Now  consider  the  <111>  situation  as  shown  in  Fig.  19.  Between  layers  n 
and  n+1,  n+2  and  n+3,  n+4  and  n+5,  etc., the  atoms  are  connected  by  three 
bonds  while  there  is  only  one  bond  between  planes  n+1  and  n+2,  n+3  and  n+4, 
n+5  and  n+6,  etc.  Based  on  this  relationship  we  assume  that  the  polishing 
removes  atoms  such  that  the  single  rather  than  triple  bonds  are  broken, 
i.e.,  all  the  plcuies  up  to  and  including  n+1  are  removed.  Thus,  the  atoms  of 
layer  n+2  will  be  pulled  toward  tne  n+3  plane.  Again,  if  the  former  atoms 
are  displaced  completely  onto  the  latter  the  bond  distance  become  ( l/fs')  a, 
corresponding  to  a  compressive  strain  of  about  6%. 

The  above  analysis  indicates  that  the  polishing  should  produced  a 
compressive  strain  and  that  the  ratio  on  the  <100>  to  <111>  surfaces  should 
be  a  factor  of  three.  Both  these  deductions  are  in  good  agreement  with  our 
experimental  results.  However,  the  absolute  strain  predicted  by  assuming 
complete  displacement  is  significantly  larger  than  our  observations.  For 
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example.  For  the  <100>  surface  18%  (predicted)  to  2-3%  (observed).  This 
indicates  that  there  may  not  be  complete  displacement  or  that  our  model  is 
too  simple  to  account  for  these  quantitative  details, 
b.  Diffusion  Model 

From  Fig.  16  and  17  as  well  as  Tables  II  and  III  we  note  for  GaAs  <100> 

and  <111>  that  d^  is  a  function  of  polishing  time  t  and  that  for  both 

situations  it  saturates  with  increasing  time.  This  time  dependence  of  d  can 

s 

be  explained  on  the  basis  of  a  diffusion  model  which  will  yield  a  diffusion 
constant,  <,  for  a  given  polishing  condition. 

We  assume  that  there  are  two  types  of  processes  occurring  during  the 
polish  procedure.  The  first  is  the  strain  is  created  at  the  surface  and 
diffuses  into  the  material  perpendicular  to  the  surface.  The  second  process 
is  the  removal  of  the  surface  layers  during  the  polish  operation.  The 
surface  strain  provides  a  driving  force,  so  that  the  strain  is  propagated 
into  the  sample.  For  simplicity,  we  assume  that  the  diffusion  occurs  in  only 
one  dimension,  that  is: 

0&a)/3t)=:K(3^Aw/9x'^)  (23) 

where  Au  is  the  shift  in  the  LO  phonon  frequency  from  its  unstrained  value, 
t  is  the  polish  time,  x'  is  the  spatial  coordinate  perpendicular  to  the 
surface  whose  origin  is  at  the  surface  when  t=0,  and  k  is  a  diffusion 
constant.  We  can  express  the  removal  of  surface  layers  during  the  polishing 
in  terms  of  the  equation, 

x'=x-t-Ct  (24) 

where  C  is  the  rate  of  surface  removal  and  x  is  the  distance  from  surface  at 
time  t.  Combining  Eqs.(15)  and  (16),  we  get  the  differential  equation: 

(3A(i)/3t)  =tc(  3^At*i/3x^)-*-C(3Aa)/3x)  (25) 

Our  boundary  conditions  are: 

Aw(x,0)=0  x>0  (26a) 
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From  experiment,  when  t>2.5  minutes,  the  strain  at  the  surface  is  almost 
independent  of  the  time  (see  Tables  II  and  III)  so  that  we  can  write: 

Auj(0,t)=  t>2.5  minutes  (26b) 

In  order  to  solve  Eq.(12)  we  make  the  addition  assumption  that: 

Au)(0,t)=  for  t>0  (26c) 

Also  from  experiment, 

limAa)(x,t)=0  (26d) 

-t 

From  the  above  boundary  conditions  the  solution  to  Eq.(12)  can  be  written 
as : 


Aw=Ato  f(x;t,d  ) 
o  '  m 


(27a) 


r 


where 

f(x;t,d  )=exp(-x/d  ){l-(2/  ix)exp(x/2d  )  \  dnexp[ -n^-x^/(  16d^^r)^ )  ] }  (27b) 


m 


m 


m 


and 


J 


m 


d  =tc/C 
m 


(27c) 


It  is  clear  that  the  mathematical  form  of  Eq.(27a)  is  different  from 
Eq.(9a),  the  expression  that  we  have  used  to  interpret  our  data.  We  will 
show  that  the  difference  between  the  two  forms  is  negligible  for  t>2.5 
minutes,  i.e.  the  shortest  time  used  in  this  work.  That  is,  we  will  demon¬ 
strate  the  Eq.(27b)  can  be  approximated  by  a  simple  exponential  (for  t>2.5 
minutes) : 

of  ^ 

f(x;t,d  )^exp[-x/d  (t)]  (28) 

ro  s 

In  the  limit  of  long  polishing  time  it  can  be  shown  from  Eq.(27b)  that 
Eq.(20)  is  rigorous.  That  is: 
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(29a) 


lim  f{x;t,d  )=exp(-x/d  ) 

and  hence 

d  =d  (00)  (29b) 

ms  '  ' 

Note  that  in  Figs.  16  and  17  as  the  time  increases  d  saturates,  in 

s 

general  agreement  with  the  above  consideration.  Hence,  from  the  experimental 

saturation  value  of  d  we  can  evaluate  d  . 

s  m 

In  order  to  plot  f(x;t,d  )  it  is  necessary  to  jcnow  not  only  d  but  also 

m  m 

the  parameter  C.  To  evaluate  this  quantity  we  polished  GaAs  <100>  and  <111> 
surfaces  with  Q.OS  pm  grit  for  30  minutes  and  measured  the  weight  loss  of 
the  samples  which  had  given  surface  areas.  From  this  procedure  we  obtain: 

C^*^*^=(6.  l+0.2)xl0  ^cm/sec 

111  -4 

C  =(5.1+0.2)xl0  cm/sec 

Using  the  value  of  given  above  and  estimating  d^=200A  from  Fig. 6, 

we  have  calculated  f(x;t,d  )  as  a  function  of  x  for  t=2.5  minutes;  these 

m 

results  are  plotted  as  the  solid  line  in  Fig.  20.  In  order  to  determine  the 

validity  of  Eq.(28),  we  have  made  a  least  square  fit  of  the  right  hand  side 

of  Eq.(28)  (a  simple  exponential  form)  to  f  (x;2.5min,200A) .  This  fit  is 

displayed  as  the  solid  line  in  Fig. 20.  The  least  square  fit  the  simple 

exponential  is  shown  by  the  dashed  line  Fig. 20.  There  is  relatively  little 

difference  between  two  curves.  This  fitting  procedure  yields 
eff 

d^  (2.5min)=110A.The  chi-square  in  this  fitting  procedure  is  about  2.1, 

indicating  that  fit  is  quite  good.  From  the  use  of  the  simple  exponential 

[Eg. (la)]  in  our  original  analysis  we  have  obtained  d  =100A  for  this  polish- 

eff  ^ 

ing  condition  (see  Table  II).  Thus  d  ^jd  and  hence  there  is  little  error 

s  s 

in  using  a  simple  exponential  [Eg. (10))  instead  of  the  more  rigorous  form  of 
Eg. (19b)  in  evaluating  the  skin  depth  of  the  polish- induced  damage.  Obvious¬ 
ly  for  t>2.5  minutes  there  is  even  a  closer  correspondence  between  f(x;t,d^) 
and  the  simple  exponential  as  indicated  by  Eg. (29a). 

Using  the  value  of  given  above  and  estimating  d^=440A  from  Fig. 17, 

we  have  calculated  f (x;2 . 5min,440A)  as  a  function  x  for  the  <111>  case.  The 
result  is  plotted  as  the  solid  line  in  Fig.  21.  As  in  the  case  of  the  <100> 
sample,  we  also  have  made  a  least  square  fit  of  the  right  side  Eg. (28)  to 
f(x;2.5min,440A) .  The  least  squcure  fit  of  the  simple  exponential  is  shown  by 
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the  dashed  line  in  Fig. 21.  There  is  also  relatively  little  difference 
between  two  curves. 

In  order  show  that  the  analysis  based  on  the  one  dimensional  diffusion 

model  is  consistent  with  the  experiment  for  different  polishing  times,  in 

©f  f 

Fig.  16  we  have  plot  d  (t)  versus  t  from  2.5  to  60  minutes  as  evaluated 

s 

by  the  above  comparison  procedure  for  the  case  of  d  =200A  (solid  line)  and 

m 

d^=175A  (dotted  line)  for  <100>  GaAs  damaged  with  O.OSpm  grit.  The  values  of 

d  (t)  were  obtained  by  a  least  square  fit  of  Eq.(lO)  to  the  set  of  data  Aw 
s 

which  are  generated  by  Eq.(27b)  with  a  particular  value  of  d  .  We  find  that 

m 

d^=175A  yields  a  better  overall  fit  to  the  experimental  results.  This  value 

is  well  within  the  error  bars  for  the  experimentally  determined  d  (t). 

eff  ® 

In  Fig.  17,  we  have  plotted  d^  (t)  versus  t  from  2.5  to  60  minutes 

for  the  case  of  d  =440A  (solid  line)  and  d  =400A  (dotted  line)  for  <111> 

m  m  -- 

©f  f 

GaAs  damaged  with  0.05ura  grit.  The  values  of  d^  (t)  were  obtained  by  a 

least  square  fit  of  Eq.(lO)  to  the  set  of  data  Aw  which  are  generated  by 

Eq.(27b)  with  a  particular  value  of  d^.  The  dotted  line  (d^=400A)  fits  the 

experiment  better  than  solid  line  (d  =440A) .  For  t<10  minutes,  the  fitting 

m 

is  not  good  as  for  the  <100>  case  (see  Fig.  16)  This  may  be  because  the 

effect  of  inhomogeneous  strain  for  <111>  surface  is  negligible,  affecting 

the  experimental  value  of  d  .  For  long  polish  times,  the  effect  of 

s 

inhomogeneous  strain  is  probably  not  a  major  factor  and  hence  there  is  a 
better  fit. 

Using  Eq.(27c)  the  diffusion  constants  for  the  <100>  and  <111>  directions 
are  evaluated  to  be: 

K:^°°=(1.07+0.05)xl0"^cm’^/sec 

K^^^=(2.04+0.21)xl0'^cm‘^/sec 

using  d  =175A  for  the  <100>  surface  and  d  =400A  for  the  <111>  case,  respec- 
tively.  The  diffusion  constant  for  the  <111>  surface  is  about  twice  as  large 
as  <100>,  indicating  that  the  propagation  of  the  strain  along  the  former 
direction  occurs  more  readily  in  relation  to  the  latter.  In  the  spring  model 
discussed  above  the  bond-breaking  for  the  <111>  surface  is  easier  than  for 
<100>.  Therefore  the  bond-breaking  and  diffusion  models  are  in  qualitative 
agreement . 
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c.  Disorder  and  Temperature  Gradient  Effect 

It  has  been  shown  that  the  influence  of  disorder  and/or  temperature 

gradients  on  the  LO  phonon  is  to  red  shift  the  peak  position  and  asymmetri- 

27 

caily  {  low-frequency  side  )  broaden  the  line  shape  .  Introducing  such  an 

effect  into  our  equations  yields  a  poorer  fit.  We  have  measured  the  RS  below 

58 

the  LO  phonon  region  and  find  little  disorder-activated  effects  .  With 
regard  to  temperature,  we  used  the  same  power  for  both  damaged  and  undamaged 
surfaces  and  so  have  no  reason  to  suspect  a  difference  in  temperature 
gradient  in  these  cases.  We  concluded  that  the  effects  of  disorder  and 
temperature  gradient  are  relatively  minor. 

d.  Free  Carrier  Effects 

There  is  also  the  question  of  whether  or  not  our  blue  shift  (of  about 
5- 10cm  is  due  to  the  introduction  of  free  carriers^*^'^^'^^  by  the  polish 

damage.  It  is  well  known  that  the  component  of  the  coupled  plasmon-LO 
phonon  mode  is  at  a  higher  frequency  than  that  of  the  uncopled  LO  peak"^  .  If 
the  scattering  volume  near  the  polished  sample  surface  contains  a  free 
carrier  concentration  of  n«1.0xl0^\m  then  L^  will  lie  about  10  cm  ^ 
higher  than  uncoupled  LO  mode  (292  cm  ^).  Thus  it  might  be  responsible  for 
the  high  frequency  bun^  in  the  damaged  spectra  of  Fig. 11(a)  and  11(b).  This 
ambiguity  can  be  resolved  by  examining  the  TO  phonon  since  this  mode  does 
not  couple  to  the  free  carriers.  We  measured  the  TO  phonon  RS  from  the  <111> 
damaged  surfaces  and  observed  a  blue  shift  from  its  unstrained  value. 
However,  because  of  the  complexities  of  the  strain-dependence  of  the  TO 
phonon  (it  is  a  doublet we  were  not  able  to  do  a  detailed  fit.  The 
observed  blue  shift  is  consistent  with  the  strain- induced  shift  of  the  LO 
mode  and  hence  we  deduce  that  free  -carriers  are  not  a  factor  in  our  analy¬ 
sis  . 

e.  Dislocations  and  Point  Defects 

In  addition  to  the  strain  due  to  the  bond-breaking  discussed  above, 

strain  due  to  the  dislocations  and  point  defects  such  as  vacancies  and 

64 

interstitials,  may  exit  .  It  is  known  that  the  dislocations  and  point 
defects  are  created  by  polish  damage  and  that  vacancies  and  interstitials 
exit  as  intrinsic  defects.  Vacancies  are  strong  centers  of  attraction, 
producing  tensile  stress  in  the  lattice,  while  interstitials  create 
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repulsion,"  producing  compressive  stress  in  the  lattice.  The  creation  of  a 
vacancy  and  interstitial  is  simultaneous,  e.g.,  there  are  equal  amounts  of 
vacancies  and  interstitials.  Therefore,  the  macroscopic  strain  due  to 
vacancies  and-  intersititials  should  cancel  each  other.  Hence  the  strain  due 
to  the  point  defects  is  not  measurable  in  this  study.  The  strain  due  to  the 
dislocations  is  a  shear  stress.  The  dislocations  may  make  a  contribution  to 
the  strain  in  the  sample.  However,  since  we  have  not  evaluated  the  disloca¬ 
tion  density  (by  a  technique  such  as  TEM)  for  the  various  polishing  condi¬ 
tions  we  cannot  comment  on  any  possible  relationship  between  our  observed 
strains  and  the  dislocations.  This  point  certainly  warrants  further  investi¬ 
gation. 

f.  Comparison  with  Previous  Investigations 

For  the  <111>  case  our  results  are  consistent  with  those  of  other 
57  58 

investigators  '  ,  including  Evens  and  Ushioda.  As  shown  in  Table  II  we 

find  that  the  strain  skin  depth  for  both  materials  is  approximately  equal  to 
the  polishing  particle  size  (for  GaAs  ds  is  somewhat  smaller  than  the 
particle  size).  Also,  for  both  semiconductors  ds  become  larger  with  increas¬ 
ing  polishing  time.  The  phonon  line  broadening  is  caused  primarily  by  the 
inhomogeneous  strain  distribution  (both  compressive  and  tensile)  in  the 
polishing  plane,  while  the  shift  of  the  RS  peak  is  due  to  the  mean  compres¬ 
sive  surface  strain. 

For  the  <100>  case  our  findings  are  somewhat  different  from  those  of 
other  investigators^^'^^'^^'^^'^^.  Other  workers  have  concluded  that  the 
strain  penetration  depth  is  on  the  order  of  the  particle  size.  However, 
they  did  not  exeunine  the  <100>  orientation  but  reported  on  other  faces. 
Also,  some  of  the  earlier  studies  were  not  "non-destructive"  but  obtained 
their  findings  by  etching  away  controlled  amounts  of  material  from  the 
surface.  We  find  that  strain  penetration  depth  for  the  <100'>  surface  is 
considerably  less  than  the  particle  size,  particularly  in  the  0.3um  and  1 
urn  cases  (see  Table  I).  Compared  with  the  <111>  surface  case,  the  strain  in 
the  <100>  surface  is  relatively  homogeneous.  The  phonon  line  broadening  is 
primarily  due  to  the  exponential  decay  of  homogeneous  strain  from  surface. 

For  <100>  GaAs  our  study  agree  with  the  <100>  GaAs  study  of  Kuwamoto 
et.  al^^ .  They  have  used  X-ray  diffraction,  RS  and  TEM  to  investigate  work 
damage  created  during  standard  processes  used  to  fabricated  GaAs  wafers. 
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Equation  (a)  W43  used  to  evaluate  the  skin  depth  while  the  concentration  of 
dislocations  was  estimated  by  cross-sectional  TEM  analysis.  The  skin  depth 
of  damage  evaluated  by  X-ray,  Raman  and  TEM  were  found  to  be  in  quantitative 
agreement . 


E.  Excitation  Wavelength  and  Pump  Chopping  Frequency  Dependence  of 

Photoreflectance  in  Hg,  Cd  Te 

1-x  X 

We  have  measured  photoreflectance  (PR)  at  77K  in  the  vicinity  of  the 
optical  feature  from  Hg^_^Cd^Te  (0.14  <  x  <  0.3)  using  infrared  (3.39  nm) 
and  deep  blue-violet  (4579A  and  4067A)  pump  beams;  both  n-  and  p-type 
materials  were  investigated^^.  The  spectral  lineshapes  were  fit  using 
existing  theories  of  electromodulation  yielding  interband  energies  and 
broadening  parameters.  In  addition,  the  dependence  of  the  PR  intensity  on 
the  modulation  frequency  (Q^)  was  studied  in  order  to  obtain  information 
regarding  surface  (interface)  traps.  The  variation  of  PR  intensity  with  Qm 
yielded  characteristic  trap  time  constants  in  the  range  1-4  milliseconds. 

In  PR  the  electric  field  that  may  occur  at  the  surface  (interface)  of  a 
semiconductor  (or  semiconductor  structure)  is  modulated  by  the 
photo- injection  of  electron-hole  pairs  by  a  secondary  (pump)  beam  chopped  at 
frequency  Q  .  The  measured  PR  spectra  are  very  similar  to  corresponding  ER 
features. 0,71,72  been  demonstrated  that  the  mechanism  of  PR  is 


Thus  PR  features  can  be  analyzed  using 
69,73-75 
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indeed  electromodulation.  ' 

existing  theories  of  electromodulation'^' ' ' "  Information  about  the 

material  and/or  structure  can  be  obtained  from  the  details  of  the  lineshape 
(energy  positions,  broadening  parameters,  etc.)  as  well  as  the  response  to 


the  modulation  (pump)  perturbation  (pump  wavelength,  modulation  frequen- 

76-78  ,  ,  ,  71,72  ^  , 

cy,  modulation  amplitude,  etc.) 


While  consideredjle  attention  has 


been  devoted  to  lineshape  features  little  work  has  been  done  in  the  latter 
76  78 

area.  '  Recently,  the  recombination  lifetime  of  electron-hole  pairs  in 

GaAs  compositional  superlattices  was  obtained  from  the  variation  of  the  PR 
77 

intensity  with  Q  .  Amirtharaj  et  al  have  recently  reported  the  first  PR 
^  76 

study  of  Hg,  Cd  Te  and  Cd,  Zn  Te.  For  the  former  material  the  3.39  um 
1-x  X  1-x  X 

line  of  a  He-Ne  laser  was  employed  as  the  pump  excitation  with  a  power 

2 

density  of  about  5  mW/cm  . 

In  this  study  we  report  the  results  of  a  PR  study  at  77K  of  Hg  Cd  Te 

X  *  X  X 

(0.14  <  X  <  0.3)  using  as  pump  excitation  sources  the  4067A  line  of  a  Kr 
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laser  4579A  line,  of  an  Ar  Laser,  or  the  3.39  iim  line  of  a  He-Ne  laser. 

No  PR  signal  could  be  obtained  for  pump  wavelengths  in  the  visible  with  an 

energy  below  the  optical  transition.  The  spectra  were  studied  in  the 

region  of  the  E  optical  structure.  The  lineshapes  were  analyzed  using 
^  7 

currently  available  theories  of  Aspnes  and  Raccah  yielding  interband 

73  75 

energies  and  broadening  parameters.  '  In  addition,  for  several  samples 
we  have  investigated  the  dependence  of  the  PR  intensity  on  pump  chopping 

frequency,  Q  .  An  analysis  of  this  variation  yields  a  characteristic  time 

^  76 

which  can  be  related  to  slow  surface  (interface)  traps, 

1.  Experimental  Details 

The  experimental  apparatus  used  to  measure  the  PR  signal,  AR/R,  was 

similar  to  that  described  in  Refs.  71  and  72.  In  addition  to  recording  the 

output  on  a  chart  recorder  the  signals  were  also  digitized  and  stored  on  a 

computer  for  analysis.  For  the  measurements  in  the  infrared  the  3.39  pm 

line  of  a  He-Ne  laser  was  employed  while  the  deep  blue-violet  experiments 

were  carried  out  using  lines  of  an  Ar  or  Kr  laser  in  the  range  4579-4067A. 

2  2 

Typical  power  densities  were  5  roW/cm  for  the  infrared  and  100  mW/cm  for 
the  deep  blue-violet.  All  measurements  were  made  at  77K  in  the  vicinity  of 
the  E^  critical  point  of  Hg^  ^Cd^Te  (~2.3-2.5  eV).  The  deep  blue-violet 
experiments  in  this  spectral  region  presented  some  experimental  problems. 
Because  of  the  proximity  of  the  pump  and  probe  wavelengths  it  was  not 
possible  to  completely  eliminate  some  of  the  strong  pump  radiation  from 
reaching  the  detector.  Thus,  the  deep  blue-violet  spectra  contained  a 
spurious  background  in  addition  to  any  real  background  that  may  be  present 
in  electromodulation  from  Hg,  Cd 

Four  different  samples  were  studied  including  bulk,  liquid  phase 

epitaxy  (LPE)  and  molecular  beam  epitaxy  (MBE)  material.  The  bulk  samples 

were  both  n-type  (n=10^^  cm  x=0.3)  and  p-type  (p=5xl0  cm  ,  x=0.3). 

The  n-type  material  was  achieved  by  annealing  for  about  30  days  Hg 

15  -3 

overpressure.  The  LPE  material  consisted  of  an  n-type  film  (n~10  cm  )  of 
thickness  lO^m  having  x=0,22  grown  on  a  CdTe  (111)  substrate.  The  MBE 
sample,  of  thickness  S^m  on  a  CdTe  (111)  substrate,  also  was  n-type 
(n=2.5xl0^^  cm  ^)  but  had  a  nominal  Cd  concentration  of  x  =  0.14.  The 
surface  of  the  p-type  material  was  polished  by  a  procedure  given  in  Ref.  82 
The  surfaces  of  the  n-type  samples  were  not  processed. 
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2.  Experiment^  Results 

^^0.78'=‘^0.22^®  (n-type) 

The  PR  spectrum  measured  from  the  LPE  Hg  ToCd  ^-Te  sample  at  77K  with 

the  3.39  urn  pump  at  Q  =  100  Hz  is  displayed  by  the  solid  line  in  Fig.  22. 

-8 

The  maximum  noise  present  in  the  spectrum  was  j^lO  .  The  signal  exhibits  a 

typical  electromodulation  lineshape.  Figure  23  displays  the  variation  in 

the  normalized  PR  intensity,  i.e.  (AR(n^)/R]/[ AR(0)/R] ,  as  a  function  of 

chopping  frequency  in  the  range  50-1000  Hz. 

The  PR  results  using  the  4067A  pump  on  this  sample  are  shown  in  Fig. 24 

by  the  dotted  curves  for  Q  =103Hz,  600Hz  and  4035Hz.  The  spurious  back- 

m 

ground  due  to  the  pump  leakage  has  already  been  subtracted.  The  lineshapes 
in  Fig.  24  are  quite  similar  to  that  in  Fig.  22.  However,  the  amplitude  of 
the  signal  is  more  than  an  order  of  magnitude  larger.  The  chopping  frequen¬ 
cy  dependence  of  the  PR  intensity  for  this  case  is  shown  in  Fig.  25  in  the 
range  40-4000  Hz. 

b.  MBE-Hg„  (n-typ») 

In  Fig.  26  is  displayed  the  experimental  PR  spectrum  (dotted  curve)  at 
77K  in  the  region  of  the  E^  optical  transition  using  4067A  pump  wavelength 
at  f=100Hz.  The  amplitude  of  the  curve  is  comparable  to  that  in  Fig.  24. 
Note  that  the  energy  of  the  E^^  optical  feature  for  this  sample  is  lower  than 
that  of  LPE  -  Hg  -.aCd  __Te  due  to  the  smaller  Cd  concentration. 

c.  Bulk  Hg^  .^Cdg  (n-type) 

The  experimental  data  at  77K  for  this  sample  in  the  region  of  the 
transition  is  plotted  in  Fig.  27  (dotted  curve)  using  the  4067A  pump  for 
f=100HZ.  The  amplitude  of  the  signal  from  this  sample  is  somewhat  smaller 
than  those  in  Figs.  24  and  26. 

d.  Bulk  Hg^  ^Cd^  ^Te  (p-type) 

For  this  sample  the  PR  spectra  at  77K  for  Q  =84Hz,  641Hz  and  3215Hz 

m 

using  the  4067A  pump  are  shown  in  Fig.  28  by  the  dotted  curves.  The  PR 
signal  at  77K  using  the  4579A  line  of  an  Ar^  laser  is  displayed  in  Fig.  29. 
We  were  not  able  to  get  a  PR  signal  from  this  sample  (or  any  of  the  other 
materials)  for  pump  wavelengths  in  the  green-blue,  i.e.  Ar  laser  lines  in 


the  range  5145A-4880A,  even  with  power  densities  (~100  mW/cm'^)  comparable  to 

the  4067A  or  4579A  lines.  Thus,  for  pump  photon  energies  below  no  PR  was 

observed.  In  Fig.  30  we  have  plotted  [AR(Q  )/AR]/[ AR(0)/R]  as  a  function  of 

m 

for  this  sample.  Note  that  the  frequency  variation  for  this  sample  is 
less  than  that  for  the  LPE  material  (see  Fig.  25). 


3.  Ana; lysis  and  Discussion 
a.  Lineshape  Analysis 

It  has  recently  been  shown  that  electric-field  modulated  reflectivity 

line  shape,  AR/R,  in  the  low-field  regime  for  Hg  Cd  Te  can  be  fit  by  a 
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generalized  derivative  functional  form  (GDFF)  given  by:  ' 

AR/R  =  E  -  4Ao.  ^ \n,E) -4AEjL^  ^  ^n  ,E )  ]  (30) 

where  E  is  the  photon  energy,  p  is  the  number  of  spectral  features  to  be  fit 
and  Cj  is  the  amplitude  of  the  feature.  The  quantity  L^^\n,E)  is  the 
lineshape  function,  with  n  depending  on  the  dimensionality  of  the  critical 
point.  The  lineshape  function  can  be  expressed  as: 


L^”\n,E)  = 


cos  [e-(M-n)<p] 


In  Eq.  (31)  0  is  a  phase  factor,  is  the  phenomenological  broadening 
parameter  associated  with  the  E^  optical  transition  and  ip  is  given  by: 

<p  =  tan"^  [r^/(E-E^)]  (32) 

The  other  parameters  in  Eq.  (30)  are  defined  in  Refs.  73  and  74. 

The  first  terra  in  Eq.  (30)  is  the  third-derivative  functional  form 

(TDFF) ■ derived  by  Aspnes  with  and  being  related  to  the  details  of  tne 

critical  point  and  applied  electric  field.  The  second  and  third  terms  have 
been  introduced  to  account  for  the  background  observed  in  the 
electromodulation  spectra  from  highly  defective  Hg^  ^Cd^Te  and  from  the  A 
and  B  faces  of  <111>  CdTe. 

The  dashed  curve  in  Fig.  22  is  a  least-squares  fit  of  Eq.  (30)  to  the 

experimental  data  (solid  line).  The  obtained  values  of  E  and  r  are  listed 

2  >■  3  ^  ^  -1 

in  Table  IV.  For  this  san^le  we  have  found  Ao  /(nfl)  =  1.28  eV  and 

-  3  -2 

AEj^/(nQ)  =  0.054  eV  .  These  quantities  are  relatively  small,  indicating 


r 


the  material  is^  fairly  free  from  defects.  However,  the  broadening  parameter 
r  is  somewhat  larger  than  the  best  values  (~70  meV)  reported  for  Hg  Cd  Te 
at  77K  of  comparable  composition. 

For  the  rneasurements  made  with  the  4067A  pump  wavelength  there  is  a 

spurious  background  due  to  stray  pump  radiation  reaching  the  detector.  This 

background  madces  it  difficult  to  distinguish  between  GDFF  and  TDFF  fits  to 

the  experimental  curves.  For  all  data  taken  under  this  condition  we  have 

found  that  both  lineshape  forms  yield  comparable  results. 

The  solid  lines  in  Figs.  22,  24,  26  and  28  are  a  least-squares  fit  of 

Eg.  (30)  to  the  experimental  data.  The  obtained  values  of  E  and  F  also 

^  ( 2 )  ^ 

are  listed  in  Table  IV.  We  do  not  report  the  coefficients  for  L  (n,E)  and 
(n,E)  since  their  interpretation  is  not  clear  due  to  the  spurious 
background  mentioned  above. 

For  the  LPE-Hg.  ToCd  material  the  values  of  both  E  and  r  are 

about  10%  different  for  the  3.39  pm  and  4067A  pump  wavelength  measurements. 

This  may  be  due  to  sample  inhomogeneity  since  no  care  was  exercised  to 

insure  that  both  measurements  were  taken  on  the  same  spot  on  the  sample. 

The  energy  of  E^  is  consistent  with  a  x  =  0.22. 

The  MBE-Hg„  o,Cd,.  ,  .Te  sample  has  E,  consistent  with  x  »  0.14  and  a 
0.86  0.14  1 

broadening  parameter  comparable  to  the  LPE  sample. 

Both  bulk-Hg^  3^©  samples  (n-type  and  p-type)  had  essentially  the 

same  value  of  E.  .  _  The  PR  results  of  Amirtharaj  et  al  at  77K  on 

*  '  76 

LPE-Hg^  -yCd^  ^Te  has  yielded  a  similar  value  for  E^.  The  quality  of  the 

n-type  sample  does  not  appear  to  be  very  good  since  it  has  =  146  meV .  Of 

all  the  samples  studied  the  bulk  p-type  sample  appears  to  be  the  best  with 

F  =  79  meV.  However,  even  this  value  is  not  as  good  as  previously  reported 

1  76  79  73 

77K  '  (or  even  300K  )  measurements. 

b.  Pump  Frequency  Dependence  of  the  PR  Intensity 

The  dependence  of  the  PR  intensity  as  a  function  of  the  chopping 
frequency  was  studied  with  a  view  towards  obtaining  information  about  the 
dynamics  of  the  surface  electric  fields.  This  variation  can  be  accounted 
for  on  the  basis  of  the  following  considerations.  The  chopped  pump  radia¬ 
tion  can  be  considered  as  a  square  wave  source.  When  light  impinges  on  the 
sample,  electron-hole  pairs  are  created.  These  charges  are  then  free  to 
fill  the  charge  traps  and  modify  the  surface  electric  field  strength. 


36 


somewhat  similar  to  the  behavior  of  a  capacitor.  We  assume  that  these 
excess  carriers  abruptly  change  the  build-in  electric  field.  When  the  light 
is  switched  off  the  trap  population  and  the  electric  field  strength  decays 
with  a  characteristic  time  t,  causing  a  restoration  of  the  original  surface 


potential.  Thus  for  chopping  frequency  f  it  can  be  shown  that  the  Fourier 

77  fil 

component  of  the  PR  intensity,  AR(Q  )/R,  can  be  written  as: 

kw  ™  /  -tl. 


where  is  the  characteristic  time  constant  of  the  i^^  trap  state. 

As  an  appropriate  means  to  extract  the  effective  average  time  constant 
r  from  the  observed  modulation  frequency  dependence  of  the  PR  intensity  we 


have  used  a  modified  version  of  Eq.  (32):  ^ 

^f2.(o)//^  4C \ -+•  ilw. T*') 


The  term  A  in  Eq.  (33)  represents  the  contribution  from  all  the  states  with 

l/t>>1000Hz,  i.e.,  very  fast  states  in  relation  to  our  frequency  range. 

The  solid  lines  in  Figs.  23,  25  and  29  represent  a  least-squares  fit  of 

Eq.  (34)  to  the  experimental  data  points.  The  obtained  values  of  x  are 

listed  in  Table  V.  All  three  measurements  yield  time  constants  on  the  order 

13  -2 

of  milliseconds.  This  is  consistent  with  the  large  number  (~10  cm  )  of 

slow  traps  with  t  >  1  millisecond  observed  at  the  HgCdTe/oxide  interface 

82 

using  capacitance-voltage  techniques.  It  should  be  noted,  however,  that 
the  LPE  Hg„  ToCd^  -.-.Te  data  presented  in  Fig.  23  obtained  using  the  3.39  urn 
pump  radiation  contains  a  large  constant  background  which  indicates  the 
presence  of  a  significant  number  of  "fast"  traps  in  addition  to  the  x  =4.4 
msec  slow  traps.  The  frequency  response  of  the  PR  intensity  from  the  same 
sample  recorded  using  the  4067A  pump  wavelength  indicated  a  dominant  trap 
level  with  x  =  2.8  msec  (see  Fig.  23).  The  apparent  differences  observed 
with  these  two  pump  wavelengths  may  be  an  indication  of  the  variation  in  the 
details  of  the  charging  and  discharging  of  the  traps  within  the  region 
penetrated  by  the  two  pump  beams.  For  instance,  traps  within  about  lOOA  of 
the  surface  and  in  the  natrive  oxixle"  layer  will' play  ■a-more  important-iole  in 
the  case  of  the  4067A  pump  due  to  its  short  penetration  depth  (~200A). 
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F.  Hg  Cd.  Te/SiO  -Photox  Interface  Properties  Studied  by  Photo-  and 

w  •  /  U 

Bias -Induced  Charging 

We  have  investigated  the  effects  of  light  illumination  in  the  wave¬ 
length  rangel|jm-2200A  (with  and  without  dc  bias)  on  the  Hg  -,Cd  _Te 

U  •  /  U  •  J 

(n-type)/SiO  -Photox  interface  at  77K  in  the 

^  83 

metal- insulator-semiconductor  device  configuration  .  Illumination  of  wave¬ 
length  X  without  bias  produced  a  flatband  shift  towards  positive 

values.  The  saturation  value  of  has  a  pronounced  peak  at  X=2800A  and 

is  negligible  for  A>5100A.  The  light-induced  charge  is  maintained  for  at 
least  8  hours  at  77K  but  leaks  off  if  the  device  is  heated  briefly  to  300K. 

We  propose  a  model  to  explain  these  results  involving  light-induced  charging 

TM 

and  discharging  of  very  slow  trap  states  in  the  SiO^-Photox  (at  or  near 
the  interface)  with  an  energy  distribution  centered  about  4.5  eV  above  the 
oxide  valence  band.  With  a  positive  gate  voltage  charging  is  extended  to 
wavelengths  of  lum  (limit  of  our  measurements),  the  resulting  charge  quickly 
leaking  off.  This  latter  effect  is  related  to  charging  of  faster  interface 
states  that  are  apparently  different  from  the  slow  states  discussed  above. 
Passivation  of  Hg.  _  Cd  Te  by  deposition  of  a  foreign  insulator  such  as 

i  X  X 

Si02“Photox  has  become  technologically  important  due  principally  to  its 

chemical  inertness  coupled  with  the  resultant  high  quality  interfaces  [86]. 

It  is  important  to  have  information  about  interface  traps  since  they  can 

limit  the  carrier  lifetimes  in  junction  devices  or,  when  charged,  alter  the 

effective  gate  bias  of  metal-semiconductor- insulator  (MIS)  devices.  The 

Hg  Cd  Te  interface  in  the  MIS  configuration  has  been  studied  in  the  fre- 
1-x  X  7  2 

quency  range  1-10  Hz  [87-90].  Work  in  the  1-10  Hz  region  on  the 

TM 

Hg  Cd  Te/SiO_-Photox  interface  [90]  revealed  a  considerable  number  of 
1-x  X  2  13-2  -1 

interface  traps  (>10  cm  -eV  ).  Illumination  effects  in  the  very  low 
frequency  domain  (-10  ^  Hz)  on  Hg  Cd  Te  MIS  devices  fabricated  with  a 

88  X 

native  oxide  indicated  a  total  trap  density  in  the  oxide  as  high  as 
13  -2 

6  X  10  cm  .  We  have  recently  reported  electroreflectaince  results  at  77K 

TM  8  X 

on  a  Hg^  ^^e/SiO^-Photox  MIS  structure 

1.  Experimental  Details 

The  structure  used  in  this  experiment  consisted  of  the  same  MIS  device 
used  in  Ref.  81.  All  light- induced  studies  were  carried  out  at  77K.  To 
provide  illumination  we  used  light  from  a  xenon  lamp  passed  through  a  0.3 
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meter  grating  monochromator  {GCA/McPherson  model  218).  In  order  to  find  the 
photon  flux  (at  \)  the  power  of  the  monochromatic  light  was  first  measured 
at  5100A  with  a  light  power  meter  (International  Light  Power  meter  model 
1L510)  and  then- -the  calibration  was  extended  to  shorter  wavelengths  using  a 
Si  photodiode  (EG&G  model  4000B)  with  a  known  spectral  response  curve.  The 

rTngi 

charging  of  the  Hg^  ^Cd^  .^Te/SiO^-Photox  interface  was  determined  oy  C-V 
measurements  at  a  frequency  of  1  kHz,  the  amplitude  of  the  ac  probe  voltage 
at  the  gate  being  50  mV. 

2 .  Experimental  Results 

a.  Light  Illumination  at  Zero  Bias  Voltage 

The  sample  (in  the  open  circuit  configuration)  was  first  illuminated  by 
light  at  \  from  the  monochromator  for  a  given  amount  of  time  (t).  The  C-V 
measurements  were  then  tciken  in  the  dark.  After  this  procedure  a  new 
illumination  period  was  started. 

From  the  C-V  measurement  we  have  established  that  if  the  sample  was 
kept  at  77K  in  the  deirk  after  the  above  illumination  procedure,  the 
light-induced  charge  remained  unchanged  for  at  least  8  hours.  This  demon¬ 
strates  that  the  trap  mechanism  is  very  slow  with  a  characteristic  frequency 
<  10  ^  sec  Thus  after  every  new  illumination  period  we  have  actually 
measured  the  cumulative  effect  of  all  preceding  illumination  periods.  After 
reached  a  value  of  20  V  (or  at  the  end  of  the  day,  whichever  came  first) 
the  sample  was  heated  to  300K  which  restored  it  to  the  initial,  uncharged 
state.  This  was  checked  by  cooling  the  sample  down  to  77K  and  taking  C-V 
measurements  after  the  30  minute  heating  procedure. 

In  Fig.  31  we  present  a  typical  set  of  C-V  curves  from  our  sample  after 
illumination  by  2800A  light  for  different  periods  of  time.  The  value  of 

in  each  measurement  was  determined  to  be  the  voltage  corresponding  to  a 

89 

capacitance  of  375  pF 

Plotted  in  Fig.  32  is  the  light-induced  flatband  shift,  as  a 

function  of  the  total  number  of  incident  photons  per  unit  area,  P(1iaj),  for 
three  different  wavelengths.  The  inset  shows  an  expanded  version  of  the 
2800A  data.  Note  that  for  all  three  wavelengths  tends  to  saturate  as  P 

increases.  The  saturation  level  for  peaks  at  A  =  2800A  while  for  A  ^ 

5100A  we  have  found  no  significant  light-induced  charging. 
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b.  Light  Illumnation  with  Gate  Bias 

Application  of  a  negative  bias  to  the  gate  produced  no  appreciable 
change  in  the  charging  picture  described  above.  However,  the  application  of 
a  positive  bias  made  it  possible  to  charge  the  system  by  utilizing  \  > 
5100A.  The  result  of  illuminating  the  structure  by  X  =  1pm  light  is  pre¬ 
sented  in  Fig.  33.  Curve  (1)  is  the  C-V  measurement  with  no  light  exposure 
and  no  bias;  this  curve  is  the  same  as  the  t=0  data  of  Fig.  31.  The  C-V 
results  after  exposure  to  1  pm  light  for  one  hour  and  a  gate  bias  of  +10  V 
is  shown  in  curve  (2).  The  light  and  the  dc  bias  were  then  shut  off  and  the 
C-V  evaluated  after  10  minutes  in  the  dark  as  displayed  in  curve  (3).  The 
results  of  Fig.  33  indicate  that  1  hour  of  illumination  under  +10  V  bias  has 
increased  from  approximately  1  V  to  4  V.  Figure  33  shows  that  the 

charge  leaks  off  relatively  quickly  since  after  only  10  minutes  of  relaxa¬ 
tion  at  77K  in  the  dark  U,.  was  reduced  to  about  2  V  from  4  V.  This  is  in 

to 

sharp  contrast  with  the  case  of  charging  by  illumination  only  and  suggests  a 
different  charging  mechanism  when  bias  is  applied. 

3.  Discussion 

a.  Charging  by  Light  Illumination  with  Zero  Bias  Voltage 

The  results  of  Figs.  31  and  32  can  be  explained  by  the  presence  of  very 

TM 

slow  electron  traps  at  or  near  the  Hg^  ^Cd^  ^Te/SiO^-Photox  interface. 
The  fact  that  the  photo- induced  flatband  shift  has  a  maximum  for  X  =  2800A 
suggests  that  two  mechanisms  are  operative;  i.e.,  both  charging  and  dis¬ 
charging.  For  2800A  <  X  <  5100A  the  light  creates  only  a  charging  effect 
while  for  X  <  2800A  both  charging  and  discharging  occur.  In  addition,  there 
is  a  distribution  of  trap  states  so  that  as  the  photon  energy  increases  (up 
to  a  certain  point)  the  number  of  traps  accessed  by  the  charging  mechanism 
also  increases,  which  explains  why  the  2800A  light  creates  a  greater 
than  say  3300A.  With  further  increase  of  the  photon  energy  the  discharging 
mechanism  also  comes  into  play  and  causes  the  2200A  light  to  be  less  effec¬ 
tive. 

The  above  mechanisms  are  illustrated  in  further  detail  in  Fig.  34.  In 

this  figure  we  display  schematically  the  interfacial  region  of  the 

TM 

Hg^  .yCd^  2Te/Si02-Photox  system.  The  energy  gaps  of  the  two  materials  are 

denoted  as  E  ,  (Hg.  _Cd-  ,Te)  and  E  .  {Si0_).  The  zero  of  energy  is  taken 
g,l  u./  U.J  z 

to  be  the  top  of  the  Si02  valence  band.  The  electron  traps  are  distributed 
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continuously  with  distribution  D^(E)  around  energy  with  full-width  at 
half-maximum  of  6E^.  The  n-type  Hg^  -j^^q  3^6  is  slightly  depleted  as 
indicated  by  the  positive  value  of  in  the  dark  (t=0  curve  in  Fig.  1). 
For  the  sake  of  clarity  the  band  gaps  of  the  Hg^  ^Cd^  (~250  meV)  and  the 
SiO^  (9  eV)  are  not  drawn  to  scale. 

As  shown  in  the  figure  two  processes  can  occur: 

(1)  Charging  -  An  electron  can  be  photo-excited  from  the  SiO^  valence 

band  into  a  trap  state,  some  of  the  resultant  holes  moving  into 
the  Hg^  ^Cd^  .^Te  region.  This  process  can  take  place  for  photon 
energies  (hw)  less  than  or  equal  to  E  _.  We  denote  the  quantum 
efficiency  of  this  process  as  Q^. 

(2)  Discharging  -  For  photon  energies  greater  than  E  _/2  not  only 

does  charging  take  place  as  discussed  above  but,  in  addition, 
electrons  in  the  trap  states  can  be  photo-excited  into  the 
conduction  band  of  the  SiO^-  Some  of  these  electrons  can  be 
transferred  into  the  conduction  band  of  the  Hg^  -jCd^  ^Te.  The 
quantum  efficiency  of  this  effect  is  denoted  as  Q2. 

Note,  that  in  order  for  process  (1)  to  cause  charging,  the  photoexcited 
hole  must  move  into  the  Hg^  3C<3q  .^Te  against  the  action  of  electric  field 
whereas  in  process  (2)  the  excited  electron  which  produces  discharging  is 
pulled  towards  the  Hg^  t^®  electric  field.  Thus  we  expect 

process  (2)  to  have  a  higher  probability  than  process  (1),  i.e. 

We  denote  as  the  number  of  states  per  unit  area  available  to  mecha¬ 
nism  (1)  but  not  mechanism  (2)  and  as  N  _  the  number  of  states  available  to 

i  /  t 

both  (1)  and  (2).  We  can  write: 

r 


N^CKcj)  = 


J  D^(E)dE 


/  D^{E)dE 
0 


Hw  <  E  J2 
9^2' 


fiio  >  E  .,/2 
9/2' 


(35a) 


(35b) 


Ni,2C^i*))  = 


jD^(E)dE 


Tiu)  <  E  _/2 

9/2 


'hw  >  E  -72 
9/2 


(36a) 


In  Fig.  34  the_dark  region  corresponds  to  while  the  shaded  region  repre¬ 
sents  N  .  In  order  to  develop  our  model  more  fully  we  define  the  follow- 

X  /  ^ 

ing  quantities: 

n  (n  )■■  =  number  of  traps  per  unit  area  filled  in  the  N,(N,  _)  part 

11^2 

of  the  trap  distribution. 

n  =  total  number  of  filled  traps  per  unit  area. 


Therefore,  for  process  (1)  we  can  write: 


dn^/dt  =  Q^(P/t)  (N^  -  n^) 

(37) 

while  for  process  (1)  and  (2): 

^"l,2^‘^^  ^  Ql(P/t)  (N^^2  ~  '^2 

The  solutions  to  Eqs.  (3)  and  (4)  can  be 

(P/t)nj_2 

written  as: 

(38) 

n^(P)  =  N^[l-exp  (-Q^P)) 

(39) 

"1,2“''  *  <'‘l,2V<V«2>l 

and 

(40) 

n[P(1iw)]  =  n^[P(fiw)]  + 

(41) 

Equations  (39),  (40)  and  (41)  show  that  n[P(hw)]  saturates  as  P-^*.  We 

denote  the  saturation  value  of  n[PCh<D)]  as  n  (liw) ,  such  that  from 

s 

Eqs.  (39),  (40)  and  (41): 


n^Chu))  =  N^(Tiu))  +  ^  22)!  (42) 

It  is  also  convenient  to  express  n[P(hu))  in  terms  of  n  (liuj)  as  follows: 

s 

n[P(hou)]=n^[l-exp(-Q^P)]  +  ^  {exp( -Q^P)-exp[ -( ^ > 

Comparison  can  thus  be  made  between  the  experimental  results  of  Fig.  32 
and  our  model.  Assuming  for  simplicity  that  the  charging  takes  place  near 
the  insulator/ semiconductor  interface,  the  photo- induced  flatband  shift, 
can  be  related  to  the  number  of  photo-filled  traps  n[  P(1i(jo)  Jby  the 
relation 

AU^[P(fiu))]  =  |e|n(Ptfia>)l/C  (44) 

where  C  is  the  capacitance  of  the  insulator  (SiO^)  layer.  Thus  can  be 

used  as  a  direct  measure  of  the  nximber  of  filled  traps,  n.  From  Eqs. 

(39) -(41)  and  (44)  the  dependence  of  AU^,  on  P(h(D)  can  be  fit  using  as 

tb 

adjustable  parameters  N. ,  N  ,,,  Q  and  Q  .  Shown  by  the  solid  lines  in 
Fig.  32  is  a  least-squares  fit  of  the  experimental  data  to  these  equations 
for  the  three  different  wavelengths. 

The  quantities  N^,  Q^,  and  obtained  by  the  above  fitting 

procedure  have  a  fairly  large  standard  deviation.  However,  the  parameter 
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n  (flu)  can  be  obtained  with  a  high  degree  of  reliability  using  Eqs.  (43)  and 
s 

(44).  Shown  in  Fig.  35  by  the  open  circles  is  n^(^uj)  determined  from 
experiment  as  a  function  of  hoj.  Note  that  there  is  a  pronounced  peak  at 
about  4.7  eV  and  that  n^(‘fiu)  =  0  for 'flu)  <  2.4  eV.. 

The  photon  energy  dependence  of  n^(hu))  can  be  accounted  for  on  the 
basis  of  the  model  in  Fig.  34  and  Eq.  (42).  From  the  above  considerations 
(least-squares  fit  in  Fig.  2)  we  find  that  =  10  and  hence  for  our  case 
Eq.  (42)  becomes 

n^Cfiuj)  =  N^('ha))  +  (1/10)  (45) 

To  calculate  and  N^^  ^  [see  Eqs.  (35)  and  (36)],  we  have  assumed  a 

Gaussian  distribution  for  D^(E)  having  E  =  4.5  eV  and  6E  =  0.5  eV.  For 

t  o  o 

E  _  we  take  9  eV  which  is  close  to  the  band  gap  of  Si0_ .  Our  value  of  E 
9,2  2  o 

corresponds  to  the  energy  of  traps  in  thermal  SiO  produced  by  P  and  In 

90  ^ 

impurities  .  Shown  by  the  closed  circles  in  Fig.  35  are  the  values  of 

n^(‘fiu))  calculated  using  the  above  procedure.  There  is  good  agreement 

between  experiment  (open  circles)  and  theory  (closed  circles). 

b.  Charging  by  light  illumination  with  positive  bias  voltage 

As  we  mentioned  before,  application  of  positive  gate  bias  extends 
photocharging  effect  towards  lower  energies  of  incident  light  (=1.2  eV  was 
the  lowest  energy  of  our  measurements).  The  accumulated  charge  leaks  off  in 
a  matter  of  minutes  -  rather  quickly  in  comparison  with  the  preceding  case. 
The  lower  light  energy  limit  of  the  effect  suggests  that  the  carrier  genera¬ 
tion  by  light  in  Hg^  7^®  substrate  is  involved.  Also,  the  faster 

charge  leak-off  points  to  a  different  nature  for  interface  states  that  hold 
the  charge.  At  the  present  time  we  are  not  able  to  propose  a  detailed  model 
explaining  this  effect. 
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TABLE 


TABLE  II4 /Values  ofi&  and  d  for  various  polishing  conditions 

on  the  <10(^  surface  of  GaM  and  InP  obtained  from  a  least 
square  fit  of  Eq8.(l)-(11)  to  the  es^rimental  data. 


Polishing 

conditions 

Material  Particle  size 

Time 

^  (D 

d 

0 

0 

s 

(um) 

(min) 

,  -l.a 
(cm  ) 

/  -Ixb 
(cm  ) 

(A)* 

0.05 

2.5 

+8. 0+0. 5 

.^3.0 

100+25 

0.05 

5.0 

+8 . 6+0 . 5 

^3.0 

120+25 

0.05 

10 

+8 . 6+0 . 5 

^3.0 

150+25 

0.05 

30 

+8 . 0+0 . 5 

43.0 

178+25 

GaAs  0.05 

60 

+9. 1+0.5 

43.0 

200+25 

0.3 

10 

+10.7+0.2 

^3.0 

235±15 

0.3 

15 

+10.5+0.5 

S3.0 

250+30 

0.3 

60 

+11.4+0.5 

4.3.0 

360+50 

1.0 

2 

+10.4+0.5 

43.0 

530+75 

0.05- 

5 

+12.6+0.5 

^3.0 

190+25 

0.05 

10 

+12.1+0.5 

43.0 

230+15 

InP  0.05 

60 

+12.1+0.5 

43.0 

300+40 

0.3 

10 

+12.3+0.3 

43.0 

385+30 

0.3 

30 

+12.5+0.5 

43.0 

440+50 

0.3 

60 

+12.5+0.5 

43.0 

533+40 

a.  Averages  over  the  value  obtained  for  the  different  laser  excitation 
wavelengths 

b.  This  remge  of  produced  no  significant  change  in  the  fit  to  the 

experimental  data. 
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TABLE  III.  Valttss  ofAw  ,Ar  ,  d  for  various  polishing  conditions  on 
the  <lllX_surface  of  GaAs  and  InP  obtain  from  a  least -square 
fit  of  Bqs.(l)-(11)  to  the  experimental  data. 


Polish  conditions 


Material 

Particle  size 

im) 

Time 

(min) 

A  u 

o 

(an“^) 

d 

s 

(A)* 

0.05 

2.5 

+1.3+0. 5 

5. 5+1.0 

100+50 

0.05 

5 

+1.4+0. 5 

6. 2+1.0 

200+100 

0.05 

10 

+1.3+0. 5 

6. 5+1.0 

340+100 

GaAs 

0.05 

30 

+1.5+0. 5 

7.0-KL.0 

400+100 

0.05 

60 

+1.5+0. 5 

7. 3+1.0 

440+100 

0.3 

10 

+1.0+0. 2 

5. 5+1.0 

1200+200 

0.3 

60 

+1.4+0 • 5 

6. 3+1.0 

4500+500 

0.05 

10 

+2. 8+1.0 

4. 8+1.0 

350+100 

InP 

0.05 

60 

+3. 3+1.0 

7. 0+1.0 

1100+200 

0.3 

10 

+3 .6+,l .  0 

6. 0+1.0 

3400+500 

0.3 

60 

+3. 6+1.0 

7 . 0+1 . 0 

4500+^00 

a.  Self-consistent  fit  to  data  for  all  five  laser  wavelengths. 
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Table  IV.  Values  of  the  energy  gc^,  B. ,  and  broadening  parameter. 

r^,  at  77K  for  several  sanples. 

Sample 

(eV) 

(meV) 

1 

f2.43"» 

96<"> 

^^®"“^0.78^'^0.22^*  <n-type)  j 

[2.40^^^ 

110<^> 

2.36<^> 

94(b) 

BuDc-Hg^  .^Cdg  (n-type) 

2.49^^^ 

146(^) 

p2.50^^^ 

79(b) 

Bulk-HgQ  .^Cdg  (p-type)  < 

|2-53^°^ 

*0 

00 

0 

/  tt  \ 

'  3.39  pm  pvinp  measurement 

4067 A  pump  measurement 

/  c ) 

'  4579A  pun^  measurement 
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Table  V.  ~  V-alues  of  the  tine  constant  t  at  77K  for  several 
Hg^_^Cd^Te  samples. 


Sample 

( 10  sec ) 


(n-type) 


Bulk-Hg^  .yCdg  2'^*  (p-type) 


L4  . 


'2.8  +  0.5 


(b) 


1.7  +  0.3 


(b) 


3.39  pm  pimp  measurement 
4067A  pump  measurement 


53 


FIGURE  CAPTIONS 


Fig.  1. 
Fig.  2. 

Fig.  3. 

Fig.  4. 

Fig.  5. 

Fig.  6. 
Fig.  7. 

Fig.  8. 

Fig.  9. 

Fig.  10. 

Fig.  11 


Raman  spectrum  of  Al^  431^0  sample  C. 

Experimental  data  (dotted  line)  for  LO^  and  LO^  for  three 
AIq  52^®^"^  samples  and  SCM  fits  for  one  parameter  (solid 

line)  and  two  parameter  (dashed  line)  cases. 

Photoluminescence  linewidth  of  several  Al-  .oln.  c-.As/InP  samples 

0 . 4o  0 . 52 

as  a  function  of  substrate  temperature,  T  . 

s 

Experimental  data  (dotted  line)  for  LO^  and  LO^  for  three 

AIq  43lriQ  g2^s/InP  samples  and  Fano  model  fit  (solid  line).  The 

values  of  the  Fano  parameters  Q  and  F  for  the  various  modes  are 

r  F 

given. 

Relationship  between  width  F  and  asymmetry  of  the  Raman  line  in 
(a)  uc-Si  and  (b)  nc-Ge,  corrected  by  instrumental  resolutiion. 
The  points  correspond  to  san^les  prepared  by  molecular  beam 
deposition. 

(a)  Electron  diffraction  pattern  and  (b)  TEM  micrograph  for  a 
Iic-Ge  film. 

(a)  Electrical  conductivity  and  (b)  Raman  measurements  for 
molecular  beam  deposition  pc-Si  as  a  function  of  annealing 
temperatures . 

(a)  Electrical  conductivity  and  (b)  Raman  measurements  for 
molecular  beam  deposition  pc-Ge  as  a  function  of  annealing 
temperatures . 

Raman  scattering  at  4579A  from  the  <100>  surface  of  n-GaAs  sample 
782-8A  for  the  configuration  x(y,z)x. 

The  intensity  ratios  1^/1^  4579A  as  a 

function  of  impurity  concentration  n^. 

(a)  Raman  spectra  of  <100>  InP  at  5145  A.  Spectrum  A  is  from  on 
undamaged  surface.  Spectra  B,  C  and  D  are  from  damaged  surfaces 
which  were  polished  with  0.05pm  particle  size  for  5,  10,  60  min., 
respectively,  (b)  Raman  spectra  of  <100>  GaAs  at  5145  A.  Spectrum 
A  is  from  on  undamaged  surface.  Spectra  B,C  and  D  are  from 
damaged  surfaces  which  were  polished  with  0.05  pm  particle  size 
for  5,  10,  60  min.,  respectively. 
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Fig.  12 


Fig.  13. 


Fig.  14. 


Fig.  15. 


Fig.  16. 


Fig.  17. 


Fig.  18. 


(a)  Raman  scattering  from  <100>  InP:  undamaged  surface  at  5145A 
(spectrum  A)  and  material  damaged  with  0.05  pm  particle  size  for  5 
min.  at  5145,  5017,  4880,  4765  and  4579  A  (spectra  B-F).  (b)  Raman 
scattering  from  <100>  GaAs:  undamaged  surface  5145  A  (spectmm  A) 
and  material  damaged  with  0.3  jim  particle  size  for  60  min.  at 
5145,  5017,  4880,  4765  and  4579  A  (spectra  B-F). 

(a)  Raman  scattering  from  <111>  InP:  undamages  surface  at  5145  A 
(spectrum  A)  and  material  damaged  with  0.05  urn  particle  size  for 
10  minutes  at  5145,  4880  and  4579  A  (spectra  B-D).  (b)  Raman 

scattering  from  <111>  GaAs:  undamaged  surface  at  5145  A  (spectrum 
A)  and  material  damaged  with  0.3pim  particle  size  for  10  minutes  at 
5145,  4880  and  4579  A  (spectra  B-D). 

(a)  For  <111>  InP  the  variation  in  LO-phonon  linewidth  r 
(deconvoluted)  and  peak  shift  Au  as  a  function  of  optical  skin 
depth,  d^,  for  the  undamaged  case  and  material  damaged  under 
various  polishing  conditions,  (b)  For  <111>  GaAs  the  variation  in 
LO-phonon  linewidth  P  (deconvoluted)  and  peak  shift  Ato  as  a 
function  of  optical  skin  depth,  d^,  for  the  undamaged  case  and 
material  damaged  under  various  polishing  conditions. 

(a)  Comparison  of  experiment  and  fit  for  <100>  InP  damaged  with 
0.3nm  particle  size  for  10  min.  at  5145  and  4880  A  excitation 
lines,  (b).  Comparison  of  experiment  and  fit  for  <100>  GaAs 
damaged  with  0.3  |.r.:  particle  size  for  10  min.  at  5145  and  4880  A 
excitation  lines. 

Damage  skin  depth  d  (t)  versus  time  for  GaAs  <100>  damaged  with 

s 

0.05  pm  grit.  The  open  circles  are  experimental  results  from 

eff 

Table  II.  The  solid  and  dotted  lines  represent  d^  (t)  versus 
time  using  d^=200A  and  d^=175A,  respectively,  for  the  least  square 
fit  comparison. 

The  damaged  skin  depth  d  (t)  versus  time  for  GaAs  <111>  damaged 

3 

with  0.05pm  grit.  The  open  circles  are  experimental  results  from 

eff 

Table  III.  The  solid  and  dotted  lines  represent  d^  (t)  versus 

time  using  d  =200A  and  d  =400A,  respectively,  for  the  least  fit 
m  m 

comparison, 

Zincblende  structure  in  the  <100>  direction.  The  distance  between 
the  layers  in  (1/4)  a,  where  a  is  the  lattice  constant. 
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Fig.  19. 


Fig.  20. 


Fig.  21. 


Fig.  22. 


Fig.  23. 


Fig.  24. 


Fig.  25. 


Fig.  26. 


Fig.  27. 


Fig.  28. 


Fig.  29. 


Zihcblende  structure  in  the  <111>  direction.  The  distance  between 
the  n  and  n+1  layers  is  3/12  a  and  the  distance  between  the  n+1 
and  n+2  layers  is  (  3/4 )a. 

Least- square  fit  of  Eq.  (la)  to  f(x;2.5  min.,  200A)  for  <100> 

GaAs.  The  solid  line  is  f (x;2.5min. ,200A)  and  dotted  line  is  the 
exponential  form. 

Least  square  fit  of  Eq.  (9)  to  f{x;2.5  min.,440A)  for  <111>  GaAs. 
The  solid  line  is  f(x;2.5min.,  440A)  and  dotted  line  is  the 
exponential  form. 

Photoreflectance  spectrum  from  the  LPE-Hg,^  ToCd.  __Te  (n-type) 
sample  (solid  line)  at  77K  in  the  vicinity  of  the  E^  transition 
using  3.39  uni  piunp.  The  dashed  line  is  a  fit  to  Eq.  (30). 
Frequency  dependence  of  the  photoreflectance  intensity  (3.39  ^im 
pump)  of  sample  LPE-Hg^  ygCd^  22"^® ’  solid  line  is  a  fit  to 

Eq.  (34). 

Photoreflectance  spectra  from  the  LPE-Hg.  __Cd-  _-Te  (n-type) 

U  »  /o  U  •  Cc, 

sample  (dotted  curves  at  77K  using  the  4067A  pump  for  modulation 
frequencies  of  103  Hz,  600  Hz  and  4035  Hz.  The  solid  lines  are 
fit  to  Eq.  ( 30) . 

Frequency  dependence  of  the  photo ref lectance  intensity  (4067A 
pump)  of  sample  LPE-Hg„  T^Cd  -.^Te.  The  solid  line  is  a  fit  to 
Eq.  (5).  . 

Photoreflectance  spectrum  from  the  MBE-Hg^  14'^®  sample 

(dotted  curve)  at  77K  using  the  4067A  pump.  The  solid  line  is  a 
fit  to  Eq.  (30). 

Photoreflectance  spectrum  from  the  bulk-Hg^  yCd^  ^Te  (n-type) 
sample  (dotted  curves)  at  77K  using  the  4067A  pump.  The  solid 
lines  eire  fits  to  Eq.  (30). 

Photoreflectance  spectra  from  the  bulk-Hg^  ^Cd^  (P"type) 
sample  (dotted  curves)  at  77K  using  the  4067A  pump  for  frequencies 
of  84  Hz,  641  Hz  and  3215  Hz.  The  solid  lines  are  fits  to 
Eq.  (30). 

Photoreflectance  spectrum  from  the  bulk  Hg^  3"^®  (P~type) 

sample  at  77K  using  the  4579A  pump.  The  solid  line  is  a  fit  to 
Eq.  (30). 
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Fig.  30. 

Fig.  31. 
Fig.  32, 

Fig.  33. 

Fig.  34, 


Fig.  35. 


Frequency  dependence  of  the  photoreflectance  intensity  (4067A 
pump)  of  sample  bulk-Hg^  -^Cd^  (p-type).  The  solid  line  is  a 
fit  to  Eq.  (34) , 

Capacitance-voltage  curves  at  77K  for  different  time  exposures  of 
2800A  light  at  zero  bias. 

Photo-induced  flatband  shift,  ,  as  a  function  of  the  total 

2 

number  of  photons  per  cm  ,  P,  for  light  of  wavelength  3300A,  2800A 
and  220QA.  The  inset  shows  an  expanded  version  of  the  2800A 
results.  The  solid  lines  are  a  least-squares  fit  to  Eqs. 

(39)-(41)  and  (44). 

Capacitance-voltage  curves  at  77K  for  (1)  no  light  exposure,  (2) 
after  exposure  to  1  Mm  light  for  one  hour  and  a  bias  of  +10  V  and 
(3)  10  minutes  after  the  light  and  bias  were  shut  off, 

TM 

Proposed  model  of  the  Hg^  ^Cd^  ^Te(n-type)/Si02-Photox 

interfacial  region  including  a  distribution  of  electron  traps  in 

the  Si0_  centered  around  E  with  full-width  at  half -maximum  of 
2  o 

6E^.  For  the  sake  of  clarity  the  energy  gaps  of  the 

Hg^  -^Cd^  2'^®(~250  meV)  and  Si02(*'9  eV)  are  not  drawn  to  scale. 

The  dark  region  represents  the  states  per  unit  area  accessed  by 
the  charging  only  process  while  the  shaded  region  indicates  the 
region  of  charging  and  dischaurging. 

Experimental  (open  circles)  and  theoretical  (closed  circles) 
values  of  n^  (hu). 
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Photoiuminescence  Linewidth 
of  AIq  4^InQ 
vs  Growth  Temperature 
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FIG.  3 
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FIG.  4 
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400  500 

AnnMiing  Temperature  (C) 
FIG.  8 
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FIG.  9 


fig.  10 
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Raman  Intensity  (arb.  units) 


Depth,  x(A) 

FIG.  20 
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GoAs  <m> 

0.05  urn  particle  size 
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FIG.  29 


fig.  30 


72 


4 


P  (IO'®/cm^) 
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FIG.  33 


SiOt 


FIG.  34 
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FIG.  35 
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